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ABSTRACT 

Aims. Photoionization models so far are unable to account for the high electron temperature rj([Om]) implied by the line intensity 
ratio [Oiii]/14363A/[Om]/i5007A in low-metallicity blue compact dwarf galaxies, casting doubts on the assumption of photoioniza- 
tion by hot stars as the dominant source of heating of the gas in these objects of large cosmological significance. 
Methods. Combinations of runs of the 1-D photoionization code NEBU are used to explore alternative models for the prototype giant 
Hn region shell IZw 18NW, with no reference to the filling factor concept and with due consideration for geometrical and stellar 
evolution constraints. 

Results. Acceptable models for IZw 18NW are obtained, which represent schematically an incomplete shell comprising radiation- 
bounded condensations embedded in a low-density matter-bounded diffuse medium. The thermal pressure contrast between gas com- 
ponents is about a factor 7. The diffuse phase can be in pressure balance with the hot superbubble fed by mechanical energy from the 
inner massive star cluster. The failure of previous modellings is ascribed to (1) the adoption of an inadequate small-scale gas density 
distribution, which proves critical when the collisional excitation of hydrogen contributes significantly to the cooling of the gas, and 
possibly (2) a too restrictive implementation of Wolf-Rayet stars in synthetic stellar cluster spectral energy distributions. A neutral gas 
component heated by soft X-rays, whose power is less than 1% of the star cluster luminosity and consistent with CHANDRA data, 
can explain the low-ionization fine-structure lines detected by SPITZER. [0/Fe] is slightly smaller in IZw 18NW than in Galactic 
Halo stars of similar metallicity and [C/0] is correlatively large. 

Conclusions. Extra heating by, e.g., dissipation of mechanical energy is not required to explain Te([Om]) in IZw 18. Important 
astrophysical developments are at stakes in the 5% uncertainty attached to [Oiii] collision strengths. 
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1. Introduction 

The optical properties of Blue Compact Dwarf (BCD) galax- 
ies are similar to those of Giant Extragalactic Hn Regions 
(GEHIIR). Their blue continuum arises from one or several 
young Massive Star Clusters (MSC), which harbour extremely 
large numbers of massive stars. 

BCDs are relatively isolated, small-sized, metal-poor galax- 
ies (Kunth & Ostlin 2000) and may be the rare 'living fossils' 
of a formerly common population. BCDs can provide invaluable 
pieces of information about the primordial abundance of helium 
(e.g., Davidson & Kinman 1985), the chemical composition of 
the Interstellar Medium (ISM, e.g., Izotov et al. 2006), the for- 
mation and evolution of massive stars, and the early evolution 
of galaxies at large redshift. Among them, IZw 18 stands out as 
one of the most oxygen-poor BCDs known (e.g., Izotov et al. 
1999) and a young galaxy candidate in the Local Universe (e.g., 
Izotov & Thuan 2004). 

The line emission of H n regions is believed to be governed 
by radiation from massive stars, but spectroscopic diagnostics 
most often indicate spatial fluctuations of the electron temper- 
ature Te (see the dimensionless parameter f^, Peimbert, 1967), 
that appear larger than those computed in usual photoionization 
models, suggesting an extra heating of the emitting gas [e.g., 
Peimbert, 1995; Luridiana et al. 1999). Until the cause(s) of this 
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failure of photoionization models can be identified, a sword of 
Damocles is hung over a basic tool of astrophysics. 

Tsamis & Pequignot (2005) showed that, in the GEHIIR 
30 Dor of the LMC, the various T^. diagnostics could be made 
compatible with one another if the ionized gas were chemically 
inhomogeneous over small spatial scales. A pure photoioniza- 
tion model could then account for the spectrum of a bright fila- 
ment of this nebula. Although this new model needs confirma- 
tion, it is in suggestive agreement with a scenario by Tenorio- 
Tagle (1996) of a recycling of supernova ejecta through a rain of 
metal-rich droplets cooling and condensing in the Galaxy halo, 
then falling back on to the Galactic disc and incorporating into 
the ISM without significant mixing until a new H ii region even- 
tually forms. If this class of photoionization models is finally ac- 
cepted, extra heating will not be required for objects like 30 Dor, 
with near Galactic metallicity. 

Another problem is encountered in low-metallicity ('low- 
Z') BCDs (Appendix A). In BCDs, available spectroscopic 
data do not provide signatures for f^'s, but a major concern 
of photoionization models is explaining the high temperature 
re([Oiii]) infered from the observed intensity ratio r([Oiii]) = 
[Oiii]/14363/([Oiii]^ 5007+4959). Thus, Stasiriska & Schaerer 
(1999, SS99) conclude that photoionization by stars fails to ex- 
plain r([Oiii]) in the GEHIIR IZwl8NW and that photoion- 
ization must be supplemented by other heating mechanisms. A 
requirement for extra heating is indirectly stated by Luridiana et 
al. (1999) for NGC 2363. 
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A possible heating mechanism is conversion of mechani- 
cal energy provided by stellar winds and supernovae, although 
a conclusion of Luridiana et al. (2001) does not invite to opti- 
mism. A limitation of this mechanism is that most of this me- 
chanical energy is likely to dissipate in hot, steadily expand- 
ing superbubbles (Martin, 1996; Tenorio-Tagle et al. 2006). It 
is doubtful that heat conduction from this coronal gas could 
induce enough localized enhancement of in the photoion- 
ized gas {e.g., Maciejewski et al. 1996), even though Slavin et 
al. (1993) suggest that turbulent mixing may favour an energy 
transfer. Martin (1997) suggests that shocks could help to ex- 
plain the trend of ionization throughout the diffuse interstellar 
gas of BCDs, but concedes that "shocks are only being invoked 
as a secondary signal in gas with very low surface brightness". 
Finally, photoelectric heating from dust is inefficient in metal- 
poor hot gas conditions (Bakes & Tielens 1994). 

Nevertheless, the conclusion of SS99 is now accepted in 
many studies of GEHIIRs. It entails so far-reaching conse- 
quences concerning the physics of galaxies at large redshifts as 
to deserve close scrutiny. If, for exemple, the difference between 
observed and computed Ts({Om\) in the model by SS99 were 
to be accounted for by artificially raising the heat input propor- 
tionally to the photoionization heating, then the total heat input 
in the emitting gas should be doubled. This problem therefore 
deals with the global energetics of the early universe. 

After reviewing previous models for IZw 18NW (Sect. |2|, 
observations and new photoionization models are described in 
Sects. [3] & El Results presented in Sect. |5] are discussed in 
Sect. |6] Concluding remarks appear in Sect. [T] Models for other 
GEHIIRs are reviewed in Appendix A. Concepts undelying the 
new photoionization models are stated in Appendices B and C. 



2. Photoionization models for I Zw 1 8 NW 

Early models are reviewed by SS99. Dufour et al. (1988) en- 
visioned a collection of small Hii regions of different excita- 
tions. Campbell (1990) proposed to enhance r([Oiii]) by colli- 
sional quenching of [Om]/I5007 in an ultra-compact structure 
(electron density A^e - 10^ cm""*). Stevenson et al. (1993) mod- 
elled a uniform sphere of radius ~ 0.4". Fairly satisfactory com- 
puted emission-line spectra were obtained, but the model H ii re- 
gions were inacceptably compact according to subsequent imag- 
ing. Firstly, IZw 18 NW is essentially an incomplete H ii region 
shell of some 5" in diameter surrounding a young MSC, which is 
not spatially coincident with the ionized gas. Secondly, both the 
highly Te-sensitive line [O m]A 4363 and the high-ionization line 
He uA 4686 are detected throughout the whole shell and beyond. 

2.1. Recent modelling attempts 

According to SS99, a 'model' is basically a uniform, matter- 
bounded spherical shell, whose only free parameter is a fill- 
ing factor e. The hydrogen density is A^h = 10^ cm"^, in- 
spired by the electron density A^e([Sii]) derived from the ob- 
served doublet intensity ratio r([Sii]) = [Sii]/t6716/[Sii]/l6731. 
The central ionizing source is a synthetic stellar cluster, which 
fits the observed continuum flux at 3327A and maximizes the 
nebular Hen emission. The inner angular radius is 1.5". The 
outer radius Voai is defined by the condition that the computed 
[Om]/i 5007/[Oii]/l 3727 ratio fits the observed one. For increas- 
ing e, the material is on average closer to the source and more 
ionized, which must be compensated for by increasing the op- 
tical depth to keep [Om]/[Oii] constant, so that the computed 



Hj6 flux increases and Heii/Hj6 decreases. For e > 0.1, the shell 
becomes radiation bounded, [Om]/[Oii] grows larger than it 
should and rout becomes less than the observed value (~ 2.5"). 
Because of these trends, SS99 discard large-e models and select 
a model with e ~ 0.01, on the basis that the computed rout, H/? 
and He ii are roughly acceptabl^H. This 'best model' presents two 
major drawbacks: (1) as stated by SS99, the computed r([Oiii]) 
is too small by a highly significant factor 1.3, and (2) [S ii] and 
[O i] are grossly underestimated. 

2.2. The [Om] line problem 

Concerning r([Oiii]), SS99 note without justification that, for 
different values of A^h, "no acceptable solution is found". It will 
become clear that this cursory statement is central for concluding 
that extra heating is required. 

In response, Viegas (2002, hereafter V02) makes the correct 
point that adopting a density less than 10^ cm""* (and e = 1) can 
help improving the computed r([0 iii]). However, in the example 
shown by V02 (A^e = 30 cm^^), r([Om]) is still 10% lovB and 
r([S ii]) is somewhat off. Moreover, not only [S n] and [Oi], but 
now [O ii] as well is strongly underpredicted, in accordance with 
the analysis of SS99. V02 then proposes that radiation-bounded 
filaments with density 10"* cm"^ are embedded in the low-density 
gas at different distances from the source. If the emission of [S ii] 
and [Oi] (together with [On]) can indeed be much increased 
in this way, this denser component entails serious difficulties. 
Firstly, since the computed [O iii]/[0 n] ratio is not very much 
less than the observed one in these filaments, a sizeable frac- 
tion of [Om] must come from them (together with [On]) and 
since, due to enhanced Hi cooling (Appendix IB.lK r([Om]) is 
now half the observed value, any composite model accounting 
simultaneously for [Oiii]/I5007 and [Oii]/l3727 will underpre- 
dict r([0 ni]) in the same manner as the uniform model of SS99. 
Secondly, since at least half the [S ii] emission should come from 
the filaments, in which r([Sii]) is again only half the observed 
value (large A^e), the composite r([S n]) will be inacceptably off. 
Thirdly, no explicit solution is exhibited and it is unclear how a 
composite model of the kind envisaged by V02 will simultane- 
ouly match all lines. From the evidences she presents, V02 is not 
founded to claim that "pure photoionization can explain IZw 18 
observations". The inconclusiveness of the alternative she pro- 
poses effectively reinforces the standpoint of SS99. 

2.3. The [O \] and [S w] line problem 

If SS99 regard the [O m] discrepancy as a highly significant fea- 
ture, they are optimistic concerning the photoionization origin 
of [Oi]/l 6300+63, underpredicted by 2dex in their model. Two 
configurations are envisaged by SS99. 

2.3.1. Extremely dense filaments? 

In a first configuration, radiation-bounded filaments of density 
10'' cm"^ are embedded in the H ii region: the density is so high 
as to severely quench most lines other than [O i] and H i: only 



' In the context, the H/3 and He ii fluxes are poor selection criteria for 
e, since the former is proportional to the (unknown) covering factor of 
the shell and the latter depends much on questionable synthetic stellar 
cluster spectra (Appendix ICt. 

^ The tolerance of 9% allowed by SS99 (adopted by V02) is probably 
too large (Sect. l3.3T >. This error bar is justified in the logics of SS99, who 
aim to demonstrate an absence of solution. 



D. Pequignot: Heating of blue compact dwarf galaxies: gas distribution and photoionization by stars in IZw 18 



3 



~ 10% of [Sn]/1 67 16+31 arises from these filaments, ensur- 
ing that r([Sii]) is not much influenced. This attempt to solve 
in anticipation the problem met by V02 (Sect. 12.21 ) raises three 
difficulties, however: (1) condensations which contrast in den- 
sity by a factor of 10"* with their surroundings and present a 
large enough covering factor (~ 10% according to SS99) as to 
intercept a significant fraction of the primary radiation would 
probably represent most of the mass; (2) since the main body of 
the model H ii region produces only one quarter of the observed 
[S ii]/16716-h31 flux, it is not clear where this doublet would be 
emittecQ and (3) this highly artificial, strictly dual density dis- 
tribution is not the schematic, first-approximation representation 
of some more complex reality, rather it is a completely essential 
feature of the model since any material at intermediate densi- 
ties would usefully emit plenty of [S n]A 6716+31, but lead to a 
totally wrong r([S ii]) as in the description by V02. 

2.3.2. Extremely distant filaments? 

The second configuration proposed by SS99 involves radiation 
bounded '[Oi] filaments' of density 10^ cm"^, located at ~ 20" 
from the source. If the spectroscopic objections of Sects. |272l & 
12.3. H are now removed since density is moderate and ionization 
is low in the filaments, new difficulties arise, notably with ge- 
ometry: (1) the filaments observed at > 10" or more from the 
NW MSC of IZw 18 have such a low surface brightness as to 
contribute negligibly to the brightness of the main shell (if they 
were projected upon it); (2) the spectrum of this weak emission 
up to ~ 15" - "Halo" of Vi'lchez & Iglesias-Paramo (1998), 
"Ha Arc" and "Loop" of Izotov et al. (2001) - shows a flux ra- 
tio [O ui]A 5007/[O ii]A 3727 of order unity, whereas this ratio is 
1/300 in the putative [Oi] filaments, suggesting that the bulk of 
the emission observed at these distances arises from a gas, whose 
density is much less than 10^ cm"-*; (3) accepting all the same 
the existence of distant [Oi] emitting regions, a very peculiar 
geometry would be required to project these regions precisely 
and uniquely upon the material of the irregular bright NW H ii 
shell to be modelled; and (4) in projection, this shell appears as 
a 1.5-2.5" 'ring', which intercepts 1/200 of a 20"-radius sphere 
(including both the front and rear sides), incommensurable with 
the covering factor ~ 1/10 assumed by SS99. 

2.4. Previous models: conclusion 

Attempts to model IZw 18 fail to explain not only r([Oiii]) but 
the [Oi] and [Sii] lines as wefl. It is diflicult to foflow SS99 
when they claim that they are "not too far from a completely 
satisfactory photoionization model" of IZw 18. The explanatory 
value of their description is so loose as to jeopardize any infer- 
ence drawn from it, including the requirement for extra heating 
in IZw 18. 

In AppendixlAl a review of models obtained for other 
GEHIIRs reveals general trends and problems, which can be 
valuably analysed using the example of IZw 18. 

3. Observations of I Zw 18 

3.1. Basic properties 

Two bright regions 5" apart, IZw 18NW and SE, correspond 
to two young MSCs associated with two distinct GEHIIRs, sur- 
rounded by a common irregular, filamentary halo of diffuse ion- 

^ SS99 prefer casting doubts on the ionization balance of S*. 



ized gas (e.g., Izotov et al. 2001), immersed in a radio Hi 21 cm 
envelope rotating around the centre of mass located in between 
the GEHIIRs (e.g., van Zee et al. 1998). Although the H i column 
density peaks in the central region, large H i structures have no 
stellar counterparts. A fainter cluster, 'Component C, deprived 
of massive stars (no prominent H ii region), appears at 22" to 
the NW of the main body. The two young MSCs, l-5Myrs 
old, are the recent manifestations of a larger starburst, which 
started some 15Myrs ago in Component C and 20Myrs ago in 
the central region (Izotov & Thuan 2004, IT04). A 20-25 Myrs 
age is consistent with the dynamics of the superbubble studied 
by Martin (1996). In a radio study, Hirashita & Hunt (2006) sug- 
gest 12-15 Myrs. IZw 18 is classified as a 'passive BCD' (e.g., 
Hirashita & Hunt 2006), that is, the MSCs themselves are rela- 
tively diffuse, the stellar formation rate (SFR) is relatively low 
(Sect. |672] i and the starburst is not instantaneous. 

That a background population 300-500 Myrs old may be the 
first generation of stars in this galaxy (Papaderos et al. 2002; 
IT04) is contested by Aloisi et al. (2007). The extended optical 
halo of IZw 18 is mostly due to ionized gas emission. Unlike 
for usual BCDs, the bulk of the stars in IZw 18 is highly con- 
centrated, suggesting perhaps a young structure (Papaderos et 
al. 2002). The distance to IZw 18, first quoted as lOMpc, has 
been revised to ~ 13Mpc (Ostlin 2000) after correcting the 
Hubble flow for the attraction of the Virgo cluster From AGE 
star magnitudes, IT04 obtain 14+1.5 Mpc. At the distance D 

- 4.0x10^'' cm (12.97 Mpc) adopted here, the diameter of the 
bright region IZw 18 NW (5") is over 300pc. From new deep 
HST photometry revealing a red giant branch and Cepheid vari- 
ables, Aloisi et al. (2007) obtain 18±2Mpc. Except for scaling, 
present results are just marginally changed if this larger distance 
is confirmed (Sect. |678T l. 

3.2. Absolute H i line fluxes and reddening 

According to Cannon et al. (2002, CSGD02), the absolute H/3 
fluxes in the 5 polygons paving the NW region and the 7 poly- 
gons paving the SE region are 4.9 and 1 .7 respectively in units of 
10 ''^ergcm"^ s '. Polygons NW D6 and SE D8 do not exactly 
belong to the main body of the H ii regions and are dismissed. 
Tenuous emission around the polygons is also neglected. 

The excess over the Case B recombination value of the ob- 
served average Ha/HyS ratios, 2.94 and 2.97 in the NW and 
SE respectively, is attributed to dust reddening by CSGD02, 
who rightly doubt the large Hi coUisional excitation obtained 
by SS99 (Appendix ID.ll l. It remains that the Balmer decrement 
is influenced by collisions and that the reddening correction 
to the observed spectrum of IZw 18 has been overestimated. 
CoUisional excitation results from a subtle anticorrelation be- 
tween Te and N(H")/N(H+) within the nebula and can only be de- 
termined from a photoionization model (Contrary to a statement 
by CSGD02, the maximum effect does not correspond to the 
hottest gas). The usually adopted recombination ratio is HafH/3 

- 2.75+0.01 (Izotov et al. 1999). It is anticipated that, according 
to present models (Sect. |5j, a better Ra/HjS is 2.83+0.02. For 
use in the present study, published dereddened intensities (also 
corrected for stellar absorption lines) have been re-reddened by 
AE(B-V) = -0.04 (in view of final results, a more nearly ac- 
curate correction could be AE(B-V) - -0.03). Then the typical 
E(B-V) for IZw 18NW shifts from 0.08 to 0.04, out of which 
the foreground Galactic contribution is about 0.02 (Schlegel et 
al. 1998). The reddening corrected H/? fluxes for the main NW 
and SE H n regions are I(Hy6) = 5.6 and 2.0 respectively in units 
of lO^'^^ergcm"^ s"'. 
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In these units, the Ha flux is 33.0 over the central 13.7x10.5" 
HST field (Hunter & Thronson 1995; CSGD02) and 42.0 over 
a 60x60" field (Dufour & Hester 1990). Adopting overall av- 
erages Ha/H/? = 2.8 and E(B-V) = 0.06, the total dereddened 
Hj6 flux for IZw 18 is 18.3. Assuming that all of the ionizing 
photon sources belong to the bright NW and SE MSCs and that 
IZw 18 is globally radiation bounded, the fraction of photons 
absorbed in the two main Hii regions is 0.41. Let Q and 2abs 
be respectively the number of photons (s"') emitted by the MSG 
and absorbed by the main shell of IZw 18NW alone. The frac- 
tion QahJQ may be smaller than 0.41 for two reasons. Firstly, 
as expansion proceeds, the shells around the starbursts become 
more 'porous' due to instabilities and the more evolved NW shell 
may be more affected. Assuming that no photon escape from the 
SE shell leads to a minimum 2abs/2 — 0.34. A more realistic 
value is probably QahJQ - 0.39±0.02, since a complete absorp- 
tion in the SE would result in a strong asymmetry of the diffuse 
halo, which is not observed. Secondly, photons may escape from 
IZw 18. This effect is probably weak, given the amount and ex- 
tension of Hi in IZw 18. The adopted nominal absorbed fraction 
for the NW shefl wifl be gabs/G = 0.37+0.03, with 0.30 a conser- 
vative lower limit, obtained for a 25-30% escape from IZw 18. 

3.3. Spectroscopic observation summary 

The optical spectrum of IZw 18 has been observed for decades 
(Sargent & Searle 1970; Skillman & Kennicutt, 1993, SK93; 
Legrand et al. 1997; Izotov et al. 1997, 1997b; Izotov & Thuan, 
1998; Vilchez & Iglesias-Paramo, 1998, VI98; Izotov et al. 
1999, ICF99; Izotov et al. 2001; Thuan & Izotov 2005, TI05; 
Izotov et al. 2006) with many instruments (hale, kpno, mmt, 
KECK, CFHT, etc), the UV spectrum with iue (Dufour et al. 1988) 
and HST (Garnett et al. 1997; Izotov & Thuan, 1999, IT99), the 
IR spectrum with spitzer (Wu et al. 2006, 2007), and the radio 
continuum with vla (Hunt et al. 2005; Gannon et al. 2005). 

The IUE aperture encompasses all of the bright regions. In 
addition to G iii]A 1909, there are indications for the presence of 
CivA 1549 and Siiii]/l 1883+92. The hst spectrum allows a di- 
rect comparison of G iii] with optical lines, but corresponds to 
such a limited area (0.86") as to raise the question of the rep- 
resentativeness of the observation for IZw 18NW as a whole. 
Nonetheless, Giii]/Hy6 is identical within 10% in available mea- 
surements, once the re-evaluation of the Ha flux within the iue 
aperture is taken into account (Dufour & Hester 1990). 

The high-resolution mid-IR spectra of IZw 18 (Wu et al. 
2007, Wu07) are secured with a 4.7x11.3" slit. Over the 
13.7x10.5" HST field, the de-redenned U/3 flux is 13.4, while the 
flux from strictly the two central Hii regions, which fill only 
part of the spitzer slit, is 7.6. The adopted Hfi flux correspond- 
ing to the mid-IR spectra is taken as 10+1, the value also used by 
Dufour & Hester (1990) for the (partial) iue aperture. Measuring 
line fluxes on the published tracings show excellent agreement 
with tabulated values, except for [S iii] 18.7//, whose flux is ten- 
tatively shifted from 2.3 to 2.8x10"'^ erg cm"^ s"' . The UV and 
mid-IR spectra are not fully specific to IZw 18NW. 

An average de-reddened emission line spectrum for 
I Zw 1 8 NW, close to the one secured by IGF99 in the optical 
range, is presented in Gol. 2 ('Obs.') of Table [3] (Line identifica- 
tions in Gol. 1 ; Gols. 3-6 are presented in Sect.|5]l. This spectrum 
differs Httle from those by Izotov & Thuan (1998) and SK93. A 
rather deep, high-resolution red spectrum is presented by SK93. 
A few weak lines are taken from a deep blue MMT spectrum by 
TI05, who however quote an [On] 3727 flux larger than in ear- 
lier studies. Absolute fluxes for H/3 and the radio continuum are 



given on top of Table [3] The 21 cm and 3.6 cm fluxes, obtained 
from Gannon et al. (2005) as a sum of 3 contours for the NW 
shell, partly originate in non-thermal processes, not considered 
here. Line intensities are relative to ¥1/3 - 1000. The intensity ra- 
tio [Nii]/i6584/.i6548 quoted by SK93 is smaller than the the- 
oretical value: this is presumably due to the presence of a broad 
Ha component (VI98). Gorrecting for the pseudo-continuum, 
the theoretical ratio is recovered and a new, smaller value is ob- 
tained for the sum of the [N ii] doublet. [O ii] 7320+30 (SK93) is 
uncertain and difficult to link to Ha. Taking into account weak 
(undetected) lines, such as [Neiv]/14724, [Feiii].i/i4702, 4734, 
4755, a continuum slightly lower than the one adopted by TI05 
leads to a moderate increase of the [Ariv] line fluxes. Lines 
[Feiv] 4906, [Feii] 5158 and [Fe vi] 5176 are seen in the tracing 
by IT05, with tentative intensities 3, 2 and 2 (H/3 = 1000) respec- 
tively. Only [Feiv] is considered in Table [3] (It is noted that the 
predicted intensities for these [Fe ii] and [Fe vi] lines will be ~ 1 ). 

The most critical (de-reddened) line ratio is r([Oiii]) = 
0.0246, the value also adopted by SS99. This is 3. 1 % larger than 
the often quoted value by SK93 (2"sHt), 0.6% smaller than the 
value by IGF99 ( 1 .5"slit) and 2.3% smaller than in the blue spec- 
trum by TI05 (2"slit). 



4. New photoionization models for IZw18NW 

Models are computed using the standard photoionization code 
NEBU (Pequignot et al. 2001) in spherical symmetry with a central 
point-like source, suited to the apparent geometry of I Zw 1 8 NW 
since the bulk of the stars of the NW MSG belongs to a cavity 
surrounded by the GEHIIR shell. Radiation-bounded filaments 
embedded in a diffuse medium are modelled. The reader is re- 
ferred to Appendices |B]&|C] for a perspective to the present ap- 
proach. Atomic data are considered in AppendixiDl 

4. 1 . Steiiar ionizing radiation 

The central source Spectral Energy Distribution (SED) is treated 
analytically, with no precise reference to existing synthetic stel- 
lar cluster SEDs (AppendixICll. No effort is done to describe 
the optical+UV continuum. The continuum flux at /13327A (de 
Mello et al. 1998) is not used to constrain the power of the MSG. 
Here, this constraint can be replaced to great advantage by the 
fraction of ionizing photons absorbed in the shell (Sect. [T2] i. 

The continuous distribution of stellar masses most often re- 
sults in an approximately exponential decrease of flux with pho- 
ton energy from 1 to 4ryd in the SED of current synthetic MSGs 
(e.g., Luridiana et al. 2003). The sum of two black bodies at 
different temperatures can mimic this shape, yet providing flex- 
ibility to study the influence of the SED. The source of ionizing 
radiation is described as the sum of a hot black body, BBl (tem- 
perature Tl> 60 kK; luminosity LI), and a cooler one, BB2 (T2 
- 40-50kK; L2). A constant scaling factor 64 (<1), reminiscent 
of the discontinuity appearing in the SED of model stars {e.g., 
Leitherer et al. 1999) and constrained by the observed intensity 
of He 11/14686, is applied to the BBl flux at > 4ryd. The ioniz- 
ing continuum depends on five free parameters. The adopted T2 
range is reminiscent of massive main sequence stars and lower 
r2's need not be considered. A sufficiently large range of Tl val- 
ues ought to be considered, as the high-energy tail of the intrinsic 
SED is influenced by quite a few WR stars, whose properties are 
either uncertain or unknowable (AppendixICb. 
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4.2. Ionized shell 

The IZw 18NW shell extends from Ri = 2.85x10^° cm to Rf = 
4.75x102" cm (1.5" and 2.5" at D = 4.0x10^5 cm). 

In 'genuine' models a smooth small-scale density distribu- 
tion is assumed (gas filling factor e unity). The gas density is 
defined by means of the following general law for a variable gas 
pressure P, given as a function of the radial optical depth, r, at 
13.6eV: 



Table 1. Constraints on model parameters" 



P(t) 



^out "t" ^in ^out 
+ 



■ tan 



Mi) 



(1) 



This law is a convenient tool to explore the effects of the den- 
sity distribution on the model predictions. P is related to the pair 
of (Te, Nh) via the ideal gas law, with derived from solving 
the statistical equilibrium equations at each step. At the first step 
of the computation (t = 0), the initial pressure is fin, while at 
the last step {t = ~ oo) the final pressure is Pout- A smooth, 
rapid transition is obtained here by adopting /c = 30 in all com- 
putations. Eq. 1 introduces three free parameters: Pin, Pout, and 
the optical depth at which the transition from inner to outer 
pressure occurs. The picture of a filament core embedded in a 
dilute medium dictates that Pin < Pout- Each filament produces a 
radial shadow, which emits much less than the material in front 
of the filament and the filament itself, since it is only subject to 
the weak, very soft, diffuse field from the rest of the nebula. The 
shadows are neglected. 

In order to represent radiation-bounded filaments embedded 
in a low density medium (Appendix lB.il ). at least two sectors are 
needed: a 'Sector 1' with » Tc (radial directions crossing a 
filament) and a 'Sector 2' with Tn, < t^. To first order, only two 
sectors are considered. Observation shows that the Hen emis- 
sion, although definitely extended, is relatively weaker in the fil- 
aments surrounding the main shell (VI98; Izotov et al. 2001). 
This deficit of He ii, unrelated to an outward decrease of the ion- 
ization parameter since Hen is a pure 'photon counting' line 
above 4ryd, suggests instead that in no radial direction is the 
main shell totally deprived of absorbing gas. With the concern 
of reaching a more significant description, the same small Tn,(3) 
= 0.05 will be attached to the remaining 'Sector 3' required to 
make up the covering factor of the source to unity in all com- 
plete models. The emission of Sector 3, a moderate contribution 
to the He ii intensity, does not impact on conclusions concerning 
the main shell and the source. 

For simplicity, in any given run, the values of the three defin- 
ing parameters of Eq. 1 are assumed to be shared by all three 
sectors. Note that and Pout act only in Sector 1. The topology 
(Appendix lB.2b of the model shell is fully determined by speci- 
fying in addition the covering factors /™'' of Sector 1 (radiation- 
bounded) and /2°'' of Sector 2 (matter-bounded), with the condi- 
tion: 



•cov 
3 



i-(./r+/r')>0' 



(2) 



and finally the optical depth Tm(2) (< Tc) of Sector 2. The full 
model shell structure depends on six free parameters. 

Adopting the same Pi and the same parameters for P(t) in the 
three sectors and assuming that the outer radius of Sector 1 is Rf 
make the computed outer radii of other sectors to be smaller than 
Rf. If, however, one would like Sector 2 to extend up to R{ and 
perhaps beyond, models should be re-run for this sector using 
now a P(r) with Pout(2) < Pin and Tc < 1 . No significant conse- 
quences on the computed spectrum result from this change, as 
the increase of radius and the decrease of density in the outer- 
most layers of Sector 2 (say, r ~ 1) have opposite effects on the 



Parameter 


Constraints 


h(H-\) 


riajr\p\ (JNo ireeaom: b(rJ-V) = U.U4) 




No ireedom K[ = z.oSxlU cm 


jrcov 


Absolute I(H/3) = 5.6x10 erg cm s 


LI, L2 


Cluster SED; / < 1.0; 2abs/G~0.37 


Ti,T2 


04 < 1.0; (log(gHc/0~-O.5) 


Sa 


He iiA 4686 


He 


He/H = 0.08; He a 5876? 


C 


Ciii]^1909 


N 


[Nii]i6584 


r\ 
\J 


[U Illj/t jUU/ 


Ne 


[Neiii]^3869 


Mg 


Mg/Ar = 10.; Mgi]i4571? 


Al 


Al/Ar= 1.; Alnai855? 


Si 


Si/Ar= 10.;Siiii]/1 1883? 


S 


S/Ar = 4.37;<[Siii]>? 


Ar 


[Arm]/! 7 135 


Fe 


[Fe III].! 4658 


One-component constant density run (A'O, A'l): 


A'h 


r([Sn]) = [Sii]i6716/[Sn]i6731 (+) 




[On]i3727 


e 


iJf = 4.75x102" cm 


One-component model with Eq. (1) (Ml): 


e= 1.0 


No freedom 




'-([Sii]) 


Tc 


[On]^3727 




iJf = 4.75xl02"cm 


Two-component model (M2, M3, M4) added freedoms: 


Jl 


' > 0; /'■'"■ = /™'' + "' < 1 


Jl 


rn,(2) < Tc; global 3-D geometry 


r„,(2) 


Fine tuning I(H/3) for given /|.^°'' 



"Question marks attached to dismissed constraints (see text). 



'local' ionization. Also, 'improving' the artificial geometry of 
Sector 3 (a thin shell at radius Pi) by assuming a lower Pin(3) or 
else a filling e « 1 would not change at all the intensity of He n, 
while the emission of other lines from this sector is negligible. 

Although three sectors are considered. Sector 3 is of no prac- 
tical consequence for the main shell and no parameter is attached 
to it. A model based on the above description will be termed a 
'two-sector model' (Sect. l53] l. 

4.3. Model parameters and constraints 

Correspondances between model parameters and constraints are 
outlined in Table [1] The parameters are interrelated and itera- 
tions are needed to converge to a solution. The weak depen- 
dance of E(B-V) on the model Balmer decrement (Sect. [T2] i is 
neglected. The SED is not fully determined by the major con- 
straint Qahs/Q. Other constraints are in the form of inequalities, 
some are half-quantitative or deal with 'plausibility' arguments. 

One emission line is selected to constrain each elemen- 
tal abundance. In Table [1] a question mark is appended to 
those lines with unreliable intensities (Table O: the intensity of 
Mg i]A 457 1 +62 is given as an upper limit as the lines are barely 
detected (TI05) and suspected to be blended with a WR feature 
(Guseva et al. 2000); detection of Aliii/l 1855 is a simple guess 
from a tracing of the hst spectrum; Siiii]/l 1883+92 is barely 
seen in the iue spectrum and only the first component of the dou- 
blet is detected in the hst spectrum (IT99). The abundances of 
Mg, Al and Si are arbitrarily linked to that of argon (Table [TJ, 
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Sum 2 BBs: 1.00 x 40000K + 1 .00 x lOOOOOK 



Sum 2 BBs: 1.00 x 40000K + 1.00 x 80000K 



Parameters 


Run" 




Model* 




of model 


NO-Nl 


Ml 


M2 


M3 


M4 




2-3 


4 


5 


6 


7 


Central source parameters 












10. 


10. 


8. 


8. 


12. 


Ll/10'" ergs-' 


3.5-1.6 


3.5 


2.0 


1.6 


1.25 


S4 


.12-.67 


0.73 


0.83 


0.93 


0.24 


T2/10-* K 


4. 


4. 


4. 


5. 


4. 


L2/10^' erg s"' 


3.5-1.6 


3.5 


2.0 


1.6 


2.5 


log(0-51. 


1.04-.71 


1.054 


0.829 


0.779 


0.734 


-log(eHe/e) 


.462-. 447 


0.445 


0.523 


0.493 


0.539 


Ionized shell parameters 










6 


.0042-.31 


1.00 


1.00 


1.00 


1.00 


P,„//t/10'cGS 


38.-8.3 


5.01 


3.40 


2.96 


2.71 






23.2 


21.7 


254 


26.8 


Tc 




5.7 


4.9 


4.0 


4.3 


fCOV 

Jl 


1.-0.43 


0.20 


0.26 


0.23 


0.29 


Jl 






0.30 


0.50 


0.60 


Tm or T^{2) 


.96-300. 


270. 


1.21 


1.46 


1.00 


Elemental abundances by number (H 


= 10^) 






C 


60-52 


45.8 


38.8 


35.3 


47.4 


N 


8.1-5.3 


3.9 


4.1 


4.0 


3.8 





198-192 


172. 


168. 


162. 


173. 


Ne 


33-30 


26.4 


25.7 


24.8 


26.9 


S (Table[B 


3.4-3.9 


5.0 


4.3 


4.3 


5.0 


Ar 


.77-.90 


1.14 


0.99 


0.98 


1.13 


Fe 


4.8-5.8 


5.8 


6.1 


6.0 


6.5 


Mean shell properties weighted by A'jXA'io,,, except for 




Q.JQ 


.20-.43 


.200 


.343 


.380 


.426 


M,JW Mo 


.15-.92 


1.02 


1.52 


1.74 


1.79 


H+/H 


.998-.98 


.957 


.961 


.963 


.948 


C/O 


.00-.017 


.041 


.038 


.037 


.052 


0+/0 


.076-. 10 


.122 


.125 


.129 


.131 


o^Vo 


.910-.85 


.784 


.793 


.790 


.773 


03+/0 


.014-.03 


.049 


.043 


.043 


.042 


re(H+)/10''K 


1.65-1.73 


1.873 


1.859 


1.896 


1.839 




1.61-1.14 


1.040 


1.013 


1.012 


1.003 


reCO+vio^K 


1.63-1.43 


1.320 


1.315 


1.309 


1.272 


re(o2+)/io*K 


1.66-1.75 


1.911 


1.915 


1.961 


1.898 




1.76-2.1 


2.576 


2.402 


2.449 


2.479 




99-18 


17.3 


11.1 


9.8 


9.4 


WeCCVcm-' 


99-8.5 


46.7 


34.0 


37.8 


41.1 


WeCO+Vcm-^ 


99-16 


62.7 


48.0 


54.0 




WeCO^+Vcm-^ 


99-18 


15.4 


10.2 


8.9 


8.5 




100-19 


10.1 


7.4 


6.3 


5.7 


f{m 


.002-.014 


.032 


.026 


.028 


.030 




.001-.018 


.022 


.024 


.023 


.020 


f2(0+) 


.0013-.02 


.021 


.024 


.024 


.023 


f2(02+) 


.002-.008 


.013 


.010 


.010 


.010 




.0010-05 


.0008 


.0006 


.0007 


.0016 



"Constant A'h. A'O: A'h = 92cm-3; A^l: A'h = 17.0 cm" 
''A'h from thermal pressure of Eq. (1). 



assuming abundance ratios close to solar (Lodders 2003). For 
simplicity, the solar S/Ar ratio is also adopted and the computed 
sulfur line intensities can be used to scale S/H according to any 
preferred criterion (Sect. 16.71 ). Hei emission lines are blended 
with strong stellar absorption lines (e.g., ICF99). He/H is set at 
0.08 by number. 

In the lower part of Table [1] are given observational con- 
straints for the structural parameters of the shell, depending on 
assumptions. In preliminary constant-density 'runs' (A^O, Nl, 
not genuine models; Sect. 15. II ), a generalization of the approach 



flu'x>*i «P (-O.™ hv/ryd) 

■ TTn.. io9(0He/Q)--0.445 


1.00 










- 


0.10 


■uii O-'J 
.- Oal»/Q - 020'-. Nl 0.67 






: Nl OoM/0 - OAS '-. (I'SC " 






hl/(ev)\ NO 0.13 




0.01 



exp (-1.00 hv/ryd) 
ioq(0H</0)--0.525 




M2 Oobi/O - ■■ 



h(/(ev) 



Sum 2 BBs: 1.00 x 50000K + 1.00 x 80000K 



Sum 2 BBs: 2.00 x 40000K + 1.00 x 120000K 



flux^^l 



- enp (-1.05 hc/ryd) 
iog(0He/Q}--0.4g3 



M3 Qobi/0 - oja 




- e«p (-1.00 hc/ryd) 
iog(QHe/Q)-- 0.539 



Fig. 1. Spectral energy distribution (SED) for I Zw 18 NW mod- 
els. Flux in ergs"' eV"' multiplied by an arbitrary constant ver- 
sus photon energy hv in eV. The SED (solid line) is the sum 
of two black bodies (dotted lines), with a discontinuity at 4 ryd 
(Sect. 14. It. Panel (a) is common to AO, Al and Ml. Panels (b), 
(c) and (d) correspond to M2, M3 and M4 respectively. Vertical 
arrows appear at 1 and 4 ryd. Straight dashed lines of the form 
exp(-£i'/!v/ryd) are exponential approximations to the SEDs over 
the range 1-4 ryd to help eyes. The exponential expressions are 
written in the panels corresponding to the SEDs [a = 0.70, 1.00, 
1 .05 and 1 .00 in panels (a) to (d), with a a measure of the EUV 
continuum hardness). Some of the SED data provided in Table|2] 
(QhJQ, discontinuity factor ^4 and 2abs/2) are recalled in the 
panels. 




Fig. 2. Gas pressure versus optical depth r for models Ml 
(dashed line), M2 (solid line) and M3 (dotted line), governed 
by Eq. 1 with parameters in Table |2] M4, close to M3, is not 
shown. For comparison, the dash-dotted lines correspond to the 
constant density runs AO and Al. 



of SS99 (Sect. 12.1b is adopted. A one-component model (Ml; 
Sect. |5.2| i shows the influence of Eq. (1). Two-component mod- 
els (M2, M3 and M4; Sect. I5.3l l generaUze Ml according to 
Sect.lO 
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5. Results 

Input and output model properties are listed in the first col- 
umn of Table |2] as: (1) five primary ionizing source parame- 
ters (Sect. 14.11 ); (2) resulting numbers of photons (s"') Q and 
gee emitted by the source above 13.6 and 24.6 eV respectively; 
(3) four (A^O, A^l, Ml) to six (M2, M3, M4) shell parame- 
ters (Sect. |4.2| i; (4) elemental abundances; (5) photon fraction 
QdbslQ absorbed in the shell; (6) mass Mgas of ionized gas in 
units of 10^ Mq, (7) mean ionic fractions of and oxygen ions 
weighted by A^e; (8) average and weighted by A^eXMon and 
average A^e weighted by A^ion for and oxygen ions. 

The model SEDs (Sect. 14. It are shown in Fig. [T] panel (a) 
is common to A^O, A^l and Ml; panels (b), (c) and (d) corre- 
spond to M2, M3 and M4 respectively. Radiation is harder and 
stronger in panel (a) (see 'hardness coefficient' a in caption to 
Fig.lTJ. The gas pressure laws Pij), drawn in Fig.|2](parameters 
in Table |2]i, illustrate the contrast between preliminary runs and 
adopted models. 

Line identifications and observed de-reddened intensities are 
provided in Cols. 1 and 2 of Table [3] Computed intensities ap- 
pear in Col. 3 and Col. 5 for Run A^O and Model M2 respectively. 
Predictions are given for some unobserved lines (intensities are 
10 times the quoted values for Hi 1215A and 2/iv). The ratios 
of computed to observed intensities, noted WO/O' and 'M2/0' 
appear in Col. 4 and Col. 6 for A^O and M2 respectively. Ideally, 
these ratios should be 1.00 for all observed lines. 

Inasmuch as the convergence is completed, at least all lines 
which were used as model constraints (Table [T]) must be ex- 
actly matched by construction (Table [3]). For the sake of eval- 
uating the models, these lines are therefore useless. Similarly, 
'redundent' lines (Hi and Hei series, etc.), which carry no as- 
trophysically significant information in the context, as well as 
unobserved lines, can be discarded. Remaining 'useful' lines are 
listed in Cols. 1-2 of Table |4] and model intensities divided by 
observed intensities are displayed in Cols. 3-8 for NO-Nl, MI- 
MA respectively. These intensities are 'predictions' in that the 
are not considered at any step of the convergence. In Table 
[Oiii]/14363A stands out as the strongest, accurately measured 
optical line. Q^hJQ is repeated in Table |4] 



5.1. Constant density runs with filling factor: NO, Nl 

NO is a preliminary run (Col. 2 of Table |2]i in which A^h is con- 
stant and Q, Que, Ri and Rf are about as in the description by 
SS99 (corrected for the larger D). The convergence process, in- 
volving O/H, Ne/H, etc. (Table [TJ, is more complete than the 
one performed by SS99, but the differences of procedures do 
not change the conclusions. If A^h is in principle derived from 
r([S ii]), the sensitivity of r([S ii]) to A^e is relatively weak at the 
low density prevailing in the shell, while the exact value adopted 
for A^H may, in this particular structure, strongly influence the 
computed spectrum. By changing coherently A^h, e and Z^"'', the 
three constraints l(Hy6), [O iii]/[0 ii] and Rf can be fulfilled along 
a sequence. A^O is extracted from this sequence by assuming, as 
in the SS99 run, a covering factor f™'' = 1. The solution is close 
to the one chosen by SS99, with A^h =92cm"^ e = 0.0042, (ra- 
dial) T ~ 1 (= T,n(2) in Table IS and r([S ii]) only -2.1% off^ the 
observed value. A^O (Cols. 3-4 of Table[3j Col. 3 of TableSi fully 
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OL. 

2.5 
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Since [O m]A 5007-1-4959 is exactly matched, the entry [O mU 4363 
in Tables|4]and|5]is the 'normalized r([Oni])', i.e., the ratio of the com- 
puted r([0 III]) to the observed r([0 iii]). 



Fig. 3. Local electron density A^e versus shell radius R in 
IZw 18 NW models: (a) A^O, (b) A^l, (c) M2; and versus optical 
depth T at 1 ryd: (d) M2. In panels (c) and (d), the vertical arrows 
indicate the boundary of the matter-bounded sector 2 of Model 
M2. Also provided are the optical thickness 'max(T)' ('> 150' 
means 'radiation bounded'), the fraction of primary photons ab- 
sorbed Q-ibs/Q, and, where relevant (A^O and A^l), the covering 
factor /™'', the filling factor e and the departure of the computed 
[S ii] doublet ratio from observation in %. 



confirms the very large problems met by SS99 with [O in]A 4363 
and [Oi]^6300 (Sect.lTITl. 

A^O also fails in that QahalQ is half the expected value. 
Decreasing Q (SED luminosity) implies to decrease A^h (for 
[Oiii]/[Oii]) and increase e (for l(Hj6)). Decreasing A^h should 
help increasing T^, thus r([Oiii]), and the high ionization lines, 
largely underestimated in A^O. Correlatively, [Oiii]/[Oii] is re- 
stored for a larger r^, which helps increasing [Oi] and other 
low-ionization lines. By further decreasing A^h and increasing 
e, whilst fine-tuning f^"^' to keep the outer shell radius Rf (and 
I(Hy6)) and ^4 (for He 11), it is possible to further increase 
until the shell eventually gets radiation bounded. The resulting 
(unique) solution is A^l (Col. 3 of Table |2l Col. 4 of Table Ell, 
which much improves upon A^O concerning [Oi] and high- 
ionization lines, while 2abs/2 = /"" = 0.43 is just a bit large. 
Because of the much lower density, A^H = 17cm"^, the ratio 
r([S 11]) is now -1-4.7% off, worse than in A^O, yet not decisively 
inacceptable. Nevertheless, the normalized r([Oiii]), enhanced 
from 0.73 to 0.82, is still very significantly too small. This fail- 
ure of A^l is illustrated in the upper panels of Figs. [3]-|5] The runs 
of A^e and with nebular radius R are shown in Figs.[3^-0H for 
A^O and Figs. [3j3-H3 for A^l. Ionic fractions of oxygen versus 
T are shown in Figs. andlSji for A^O and A^l respectively. In 
A^O, Te is everywhere above 1. 55x10"* K. In A^l, the inner Te is 
3 800K higher than in A^O, but is abundant up to t = 15, 
where is below 1.40xlO'*K and the average re([Oiii]) is not 
much increased. 

This generalization shows that no solution with constant A^h 
exists even for the rather hard SED adopted by SS99. In an ex- 
treme variant of A^l, the SED is just one 10^ K black body (con- 
verged L = 2.7x10'*' erg s-i, Sa = 0.41, log(eHe/e) = -0.27), 
but the normalized r([Oiii]) = 0.87 is still too small, despite the 
unrealistically hard SED. 
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Fig. 4. As in Fig. |3]for local electron temperature TJIQ'^ K. The 
maximum Tg and average Tg's weighted by A^eXN(H^) and by 
A^eXN(0^^) are noted along each curve. In panel (c), <re> is 
given for both the radiation bounded sector (<sect.l>) and the 
full model M2 (<mod.>). 
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Fig. 5. Local fractional concentrations of (dotted line), 
(dot-dashed line), 0^+ (solid line), (dashed Une) and 
(dotted line again) versus t in IZw 18NW models: (a) A^O, (b) 
A^l, (d) M2; and versus R: (c) M2. Vertical arrows as in Fig. [3] 



5.2. A one-sector photoionization modei: Ml 

Model Ml includes the same primary source as Run A^O and 
again only one sector, but with A^h controlled by Eq. 1 (Fig.|2]i. 
Parameters appear in Col. 4 of Table |2] and predictions in Col. 5 
of Table H The litigious lines [Ora]/14363 and [Oi]/i6300 are 
very much improved compared to Run A^O (and even A^l), as are 
[Sii] and [Oiv]. 

The decisive merit of Model Ml is to demonstrate that, with 
no extra free parameter, no change of shell size and no significant 
change of source SED, the 'r([Om]) problem' met in NO can 
be essentially solved just by considering radiation-bounded fila- 
ments embedded in a lower density (higher ionization) medium 
instead of a clumped shell at constant density. While the nor- 
malized r([Oiii]) is 0.96 ([Oiii] in Table |4]i, no dense or distant 
clumps of the kind postulated by SS99 (Sect. |2]i are needed to 
account for low-ionization lines. The pressure contrast is < 5. 



Qdbs/Q is again too small, but Q cannot decrease because 
64 is close to unity (Table |2]l. In A^O, 64 was (perhaps anoma- 
lously) small, enabling a shift from A^O to A^l. Hardening the 
already hard primary (a = 0.7 due to large Tl; Fig. [TJ would 
enhance [Ariv], predicted too strong. Also, /™'' is only 0.20, re- 
sulting in an artificial cigar-like radial distribution, in which the 
low-density gas exactly shields the denser filaments from direct 
primary radiation. 

Obviously, the Umits of the one-sector model are being 
reached. A matter-bounded sector is to be added for the sake 
of a larger absorbed fraction of photons in the shell, but not prin- 
cipally to improve the akeady quite satisfactory intensities of 
[Oiii]/! 4363 and [Oi]/I6300. 

5.3. Two-sector pliotoionization modeis: Ml, M3, MA 

The enhancement of Hen and [Oiii] related to the matter- 
bounded sector must be balanced by a weaker/softer SED. 
Models illustrate the influence of the SED (Fig.[T]). 

Given a SED and both covering factors, then t^{2), and 
64 can be fine-tuned to account for I(Hy6), [O in]/[0 11] and He n. 
Iterations along the same lines as for Model Ml eventually lead 
to a model, provided that the limits on parameters are respected 
(Table [TJ. A two-parameter model sequence can be attached to 
any SED by considering several pairs (/™'', /j'" ), but little free- 
dom is attached to /J^"' , as Sector 1 is where ~75% of H/? and 
most of [O n] come from. Then /J^"'' must be on the order of, or 
moderately larger than the /'"'' of the one-sector model, say, in 
the range 0.2-0.3. Also /™'' cannot be small since one-sector 
models are rejected (Sect. |5.2| l and /["'' -H /""' must be kept sig- 
nificantly less than unity to enable He 11 excitation beyond the 
shell (Sect.gji, implying f^"" ~ 0.3-0.6. 

5.3.1 . Model M2 and variants 

Model M2 (Col. 5 of Table |2]l is the first 'complete' model. 
Gabs/G is at the low end of the nominal interval. Ouput of line 
intensities (Cols. 5-6 of Table [3j Col. 6 of Table |4| is to be con- 
trasted to the A^O output. Runs of A^e and with R are shown in 
Figs. |3};-@j;. The sharp 'spike' of the A^e curve shows how thin 
a radiation-bounded filament is compared to the shell. Runs of 
A^e and are best seen in plots versus r (Figs.[3}l-H}l)- In plots 
for M2, vertical arrows mark the outer boundary of Sector 2. Ion 
fractions 0"^/0 versus R and t are shown in Figs. |5]:-]5}l. 

Sufficient ionization is maintained in M2 owing to the lower 
average density, which also helps increasing in the high- 
ionization layers, despite the significantly softer radiation field 
(smaller QhJQ and larger a. Fig. [T]5): the inner T^. is now 
2.5x10"^ K. In Fig.Bl, the jumps of Te at r ~ 0.1 and -4.7 coiTe- 
spond to the boundaries of the He^^ shell (fairly well traced by 
O^*^ in Fig.|5}l) and the filament respectively. Comparing Fig.|5}l 
to Fig. [Sji, the ionic fractions are qualitatively similar in Model 
M2 and Run A^l, but the transition from O^^ to O^ is sharper 
and occurs at a smaller optical depth in M2. 

Average properties and abundances of Model M2 are quite 
similar to those of Ml and the predicted [Oiii]/l4363 intensity 
is again slighly weak, although the score of M2 is significantly 
better for [Ariv] and [S iv] (Col. 6 of Tablegli. 

Variants to Model M2 can be obtained by changing /["'' and 
Z™^' within limits, while retaining source parameters (except for 
minute fine-tuning of ^4). In Col. 1 of Tableware listed 7 shell 
parameters and 6 lines extracted from Table |4] M2 (Col. 2 of 
Table 15) is compared to re-converged models M2b (Col. 3) and 
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M2c (Col. 4). Increasing Z™'' from a small to a large value, with 
/™'' left unchanged, structure parameters are not much changed 
except for a decrease of t^{2) and a small decrease of O/H due to 
the larger weight of the hot high-ionization zone. Accordingly, 
[Ariv] is increased, but [Oiii]/l4363 is increased by only 1%. 
Decreasing /['"' from 0.26 to 0.22, H/3 is recovered by increas- 
ing Tn](2) and [On]/i3727 by decreasing Tc, with the conse- 
quence that Pin must decrease, thus increases and O/H de- 
creases. Finally the A 4363 intensity increases up to the observed 
value, [Ariv]i4740 and [Siv]^10.5jU increase and [Oi]/I6300 
decreases. 



5.3.2. Models M3 and M4 

In Model M3, 72/10"* K is enhanced from 4 to 5. gabs/fi is larger 
than in M2 due to lower luminosity. The larger T2 increases the 
average energy of photons absorbed in the O^^ region and the in- 
tensity of [Oiii]/l4363 is slightly larger. In the selected example, 
[Om] is again exactly matched as in M2c, but for more 'stan- 
dard' f""' and Z™'' (Col. 6 of Table|2]i. Line intensity predictions 
are slightly improved (Col. 7 of Table|4] versus Col. 4 of Table|5]). 

In Model M4, T2 is again as in M2, while Tl/IO'^K is en- 
hanced from 8 to 12 and LI/L2 is halved. Q^hs/Q is close to its 
allowed maximum due to radiation hardening, which also leads 
to (^4 « 1 (Col. 7 of Table|2]). The large flux just below 4ryd en- 
hances simultaneously the high and low ionization lines, but the 
heating of the O^^ region is lesser and the large Tl (positive cur- 
vature of the SED) does not favour a large r([Oiii]) (Col. 8 of 
Tablegl. From Table|5l variants M4b (Col. 5) and M4c (Col. 7) 
of M4 (Col. 6) fail to enhance [Om]/l4363 up to the observed 
value. Increasing the source luminosity by 20%, thus decreasing 
QdbJQ from 0.43 to 0.36, has strictly no effect on the predicted 
[O in] after convergence. 



6. Discussion 

Irrespective of the 'technical' demand raised by 2abs/2 in 
Sect. 15.21 a two-sector model is the minimum complexity of any 
shell topology (Sect. |4~2] i. The two degrees of freedom attached 
to the matter-bounded sector are inescapable. 

Model MA (and variants) appears slightly less success- 
ful than M2 and M3 concerning [Om]/l4363 and the high- 
ionization lines ([Ariv]/t 4740, [S iv]/l 10.5yu). MA has a possibly 
less likely SED and presents the largest Pout/Pin- The discussion 
focuses on Models M2 and M3, with M2 the 'standard' from 
which variants are built. 



6.1. Spectral energy distribution 

Accounting for a QnhJQ larger than, say, 1/3 turns out to be de- 
manding. Selected models correspond to nearly maximum pos- 
sible values for each SED. Acceptable 2abs/G can indeed be 
obtained, but the latitude on the SED and power of the ion- 
izing source is narrow. The uncertainties in evolutionary syn- 
thetic cluster models (AppendixO and in the evolutionnary 
status of IZwl8NW itself are sufficient to provisionally ac- 
cept the 'empirical' SED corresponding to preferred models M2 
or M3 (Fig. [T]i as plausible. The 'predicted' typical trend is 
Lv{hv- 1 ^4ryd) oc exp(-/iv/ryd). 



6.2. ionized gas distribution 

The model is most specific in that emission lines partly arise 
from a low density gas, while the largest A^e is A^e([S n]). A den- 
sity < 10 cm"^ (Table |2]i appears very low by current standards 
of photoionization models for BCDs (Sect. lA.lb . Nevertheless, 
the superbubble model of Martin (1996) is consistent with a cur- 
rent SFR = 0.02 Moyr ' for the whole NW-hSE complex. The 
two best estimates in the compilation by Wu07 are 0.03 and 
0.02 Mq yr"' . Adopting half the Martin (1996)'s rate and a wind 
injection radius of 0.1 kpc (1.5" at 13Mpc) for the NW cluster 
alone, expression (10) in Veilleux et al. (2005) suggests an in- 
ner pressure of the coronal gas P/k ~ 3x10^ Kcm"^, hence an 
ambient ISM number density ~ 12 cm"^ in the inner region {Ts~ 
2.5xlO'*K; Sect. IS, or else A^h - 12/2.3 = 5cm-^ for a pho- 
toionized gas in pressure balance with the coronal phase which 
presumably permeats the shell. Although e is unity in models, 
this phase can fill in the volume corresponding to Sector 3. 

6.3. Opticai and UV iines 

[Oiii] lines are discussed in Sects. |6^ & I6. 81 Civ and Sim] are 
accounted for within uncertainties. Other UV lines are elusive 
(Col. 6 of Table EJ. 

Computed fluxes for 8 optical lines are within 20% of ob- 
servation (Table |4]l, which is satisfactory considering the weak- 
ness of some of the lines. The 10-15% discrepancy on the ra- 
tio [S III] 6312/9531 does not challenge the model itself, given 
the various uncertainties. The ~ 20% underestimation of [S ii] 
should be considered with respect to [Siii]. The line /1 9531 is 
matched, but the far-red flux may be less reliable, and /1 63 12 
departs from observation about in the same way as [Sii], but 
/1 63 12 is a weak line. The exact status of [S ii] is undecided. 

The A^e-sensitive intensity ratio [Feiii]/14986//14658 is 
somewhat small in A^O, large in A^l and more nearly correct in 
Mi models. Would [S ii] be emitted in a high-A^e gas compo- 
nent as suggested by SS99 and V02, then [Fe iii] 4986, roughly 
co-extensive with [On], would be undetectable. The weak line 
[Feiv]4906, co-extensive with [Oiii], confirms the ion distri- 
bution of the Mi models and the iron abundance, although the 
agreement with observation is partly fortuitous. [Fev]4227 is 
overpredicted by a factor ~3, but the observed intensity is very 
uncertain and could be 2-3 times stronger than the quoted value, 
as judged from published tracing (TI05). Also, only one com- 
putation of collision strengths has ever been done for the opti- 
cal lines of the difficult [Fev] ion (Appendix |D.3t . Finally, the 
ionic fraction Fe^^/Fe, less than 5%, is subject to ionization 
balance inaccuracy. The predicted intensity ~ I(/l4227A)/4 of 
[Fe v] 407 1 enhances the computed flux of [S ii] 4068 up to the 
observed value (Tables[3]&|4]l. 

The observed He i line intensities are inconsistent (Table O, 
due to stellar lines (Sect. |43T l. [Ariv]4711, blended with 
He 1 47 13, is therefore useless. The weak [Ariv]4740 tends to 
be overestimated by ~ 50% in the preferred models. Trial calcu- 
lations show that, adopting a recombination coefficient 12 times 
the radiative one (instead of 8 times, Appendix |D.3t and dividing 
Ar/H by 1.13, [Arm] and [Ariv] would be matched in M2. 

6.4. Infrared fine-structure lines 
6.4.1. [Neiii] and [Sin] 

The reliably observed IR lines with optical counterparts, 
[Neiii] 15.5ju and [Siii] 18.7//, are very well matched, confirm- 



10 



D. Pequignot: Heating of blue compact dwarf galaxies: gas distribution and photoionization by stars in IZw 18 



ing the scaling adopted for the spitzer fluxes and the model tem- 
peratures. The Mi models, globally hotter, are more successful 
than the Ni runs. No in excess of the one of the adopted con- 
figuration (Fig.|8}l) is required. 

The predicted intensity of [S m] 33.5/i is only 40% of the ob- 
served value. Since the theoretical ratio of the [S iii] IR lines is 
insensitive to conditions in IZw 18, looking for alternative mod- 
els is hopeless. The collision strengths D. for the [S iii] lines may 
not be of ultimate accuracy, as the results of Tayal & Gupta 
(1999) and Galavis et al. (1995) differ, but the more recent Q's 
are likely more accurate. Also, the predicted [S iii] 33.5yU is even 
worse using older data. Since Wu07 cast doubts on the accu- 
racy of the flux calibration at the end of the spitzer spectrum, it 
is assumed that the [S iii] atomic data are accurate and that the 
observed fluxes around A 34yU should be divided by 2.3. 

6.4.2. [Oiv] 

If the drift of flux calibration at A 34// (Sect. |6.4TTl ) smoothly van- 
ishes towards shorter wavelengths, the [O iv] 26// flux may still 
be overestimated. Conversely, the spitzer field of view encom- 
passes IZw 18 SE, which emits little Hen, leading to underesti- 
mate [O iv]/Hy6 in I Zw 1 8 NW. Since these effects act in opposite 
directions, the original [O iv]/Hy6 is adopted for I Zw 1 8 NW. 

As shown in Fig. |5] O-*^ and O^^ coexist in the He^^ zone. 
Q3+/Q2+ ^jjj therefore [Oiv]^25.88///Heni4686 as well are 
sensitive to A^e- In NO, Hen is matched and [Oiv] is strongly 
underpredicted (Table |4]l. The predicted [Oiv] flux improves 
in the conditions of A^l and even by-pass observation in Ml, 
whose ionizing flux is however too large (Sect. lSSt . In the stan- 
dard Model M2, the predicted [O iv] exactly matches observa- 
tion after adding the blended line [Fe v] 25.91//, whose com- 
puted flux is ~ 2% of [Oiv]. Since relevant atomic data are re- 
liable, [Oiv] 25. 88fi indicates that Ne must be on the order of 
lOcnT^ in the Hen emitting region of the IZw 18NW shell. 
The model density results from general assumptions (photoion- 
ization by stars, shell geometry, e - 1, Eq. 1, etc.) and a require- 
ment to match a few basic line intensities with no reference to 
high-ionization lines, but Hen. The computed [Oiv] intensity 
is a true prediction, especially as the models were essentially 
worked out prior to IR observations: the spectrum presented by 
Wu et al. (2006) showed the predicted [O iv] line, finally noted 
by Wu07. 

6.4.3. [S IV] 

The predicted [Siv] 10.5// flux is twice the observed one. The 
collision strengths obtained by Tayal (2000) and Saraph & 
Storey (1999) for this line are in good agreement. The average 
fractional concentration of S"*^, ~ 1/3, is stable in different mod- 
els because sulfur is mostly distributed among the three ions 
g2+_|g4+ Displacing the ionization balance by changing, e.g., 
the gas density tends to make either S^^ or S**^ to migrate to S^^. 
Only in the unsatisfactory run A^O is [S iv] accounted for. 

Two-sector, constant-pressure 'models' allowing e<l and 
using the SED of Model M2 were run with the conditions 
QabJQ > 0.3 and O/H < 1.7x10 ^ In these trials, the [O iii] 5007 
and r([Sii]) constraints (Table [TJ are relaxed and the observed 
[Siv]/[Siii] ratio is exactly matched by playing with A^e (gas 
pressure) and e. Despite ample freedom and because of the 
higher A^e ~ 25 cm"-', the computed [O iv] flux is at most 60% of 
the observed one. Thus, forgetting other difficulties, the sugges- 
tion is that the excess [S iv] flux can only be cured at the expense 




10 100 1000 



Fig. 6. As in Fig. [TJ) for the variants of M2 comprising a hot 
black body emitting at soft X-ray energies: Model M2x (solid 
line) and M2x2 (dash-dotted line). The unit flux on the ordinate 
is now lO'*' erg s ' eV ' . The dash-triple-dotted line is one pos- 
sible fit to the CHANDRA data for I Zw 1 8 NW. 



of [Oiv]. A broader exploration of the SED (discontinuities?) 
and the gas distribution could be undertaken. 

The [S IV] 10.5// flux pubhshed by Wu et al. (2006) was 25% 
larger than according to Wu07. The new value should be pre- 
ferred, but this difference is at least indicative of possible uncer- 
tainties. The ratio [S iv]/HjS may also be intrinsically larger in 
IZw 18 NW than in IZw 18 SE. 

The theoretical ionization balance of some ions of sulfur 
(and argon) is subject to uncertainties (Appendix lD.3l l. The ob- 
served [S iv]/[S III] ratio can be recovered in M2 if the S^^ re- 
combination coefficient is multiplied by a factor 2.3, which is 
perhaps too large a correction: then, both [S iii] and [S iv] are 
matched if S/H is divided by 1.33, with the caveat that the pre- 
dicted [S ii] intensities are divided by 1.3. A combination of ob- 
servational and theoretical effects just listed could alleviate the 
' [S iv] problem' . 

6.5. Low ionization fine-structure lines 

Although the error bars of order 10% quoted by Wu07 may 
not include all sources of uncertainties, both [Nen] 12.78// and 
[Fe ii] 25.98// are detected in high-resolution mode. [Si ii] 34.80// 
is strong, even though the flux quoted by Wu07 (Table |4| may 
be too large (Sect. 16.4. 11 1. Usually, the bulks of [Sin] 35// and 
[Fe ii] 26// arise from a Photon-Dominated Region (PDR), at the 
warm H i interface between an ionization front and a molecular 
cloud (e.g., Kaufman et al. 2006). Schematically in a PDR, the 
photo-electric heating by UV radiation on dust grains (and other 
molecular processes) is balanced by fine-structure (and molec- 
ular) line emission. The small reddening intrinsic to IZw 18 
(Sect. 13. 2t and the Targe' gaseous iron content (Sect. 16.71 ) imply 
that little dust is available. Molecules and PAHs are not detected 
in IZw 18 (Vidal-Madjar et al. 2000; Leroy et al. 2007; Wu07). 
The classical PDR concept may therefore not apply to IZw 18, 
raising the question of the origin of [Si ii] 35// and [Fe ii] 26//, 
both underpredicted by factors 5-10 in the models (Table |4|i. 

A way to produce a 'pseudo-PDR' is X-ray heating. Two 
new variants of M2 are considered, in which a hot black body 
representing a soft X-ray emission from IZw 18NW is added 
to the original SED. The adopted temperature is Tx - 2x10^ K 
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Fig. 7. Model M2x for IZw 18NW. (a) Electron density A^e 
(solid line) and total hydrogen density N(H) (dotted line) ver- 
sus optical depth r at 1 ryd: N(H) levels out at 300 cm"-'. (b) 
Electron temperature TJIO'^ K: the maximum value and average 
values weighted by A^eXN(H+) and by A^eXN(0^^) are noted for 
the full 2-sector model (< mod. >). (c) Ratio of A^e to atomic hy- 
drogen density N(H"). In panels (fl)-(c), the dashed lines corre- 
spond to M2x2- (d) Local fractional concentrations of 0° (dotted 
hne), 0+ (dot-dashed line), 0^+ (solid Hne), 0''+ (dashed Hne), 
O"*"^ (dotted line again) and O^^ (solid line again): the thinner 
dash-dotted line (t > 200) corresponds to in M2x2 and the 
thinner solid line (split into Mix and Mlxi for t > 200) to Ne+. 



and the luminosity Lx = 4 and 8x10"'^ ergs"' for variants M2x 
and M2x2 respectively (Fig. |6|l. The actual X-ray luminosity of 
IZw 18, ~ 1.6x10^^ erg s"' in the 0.5-10 kev range of chandra, 
mainly arises from the centre of IZw 18NW and is consistent 
with a power law of slope -1 (Thuan et al. 2004), drawn in 
Fig. |6] The SED for M2x is a relatively high, yet plausible ex- 
trapolation of the CHANDRA data. M2x2 is considered for com- 
parison purpose. A^H, still governed by Eq. (1), is leveled out at 
300 cm-^ 

The run of physical conditions with t is shown in Fig.|2]for 
M2x. While <re> is increased by only ~ lOOK in the Hii re- 
gion, a warm low-ionization layer develops beyond the ioniza- 
tion front. Computation is stopped at = 100 K (r ~ 7x10"*, 
compared to 250 in M2). In M2x2 (dashed lines in Figs.|7^-[2};), 
the new layer is hotter and more ionized (final t ~ 1 .2x10^). The 
geometrical thickness of the H i layer is 2 1 % and 35 % of the H ii 
shell in M2x and M2x2 respectively. In the H i zone, O^/O does 
not exceed ~ lO^^*, in marked contrast with Ne^/Ne, overplotted 
as a thin solid line in Fig.|7}l. Lines [Ne ii] 12.8// and [Am] 7.0// 
are usually discarded in PDR models on the basis that the ion- 
ization limits of Ne" and Ar" exceed 1 ryd. Here, owing to the 
scarcity of free electrons and the lack of charge exchange with 
H", photoionization by soft X-rays can keep 1-10% of these el- 
ements ionized. 

In Hi regions, cooling is due to inelastic collisions with H". 
Reliable collisional rates exist for the main coolents [Cii] 157yU 
(Barinovs et al. 2005) and [Oi] 63// (Abrahamsson et al. 2007; 
several processes need be considered for [Oi]: see Chambaud 
et al. 1980; Pequignot 1990) and for [Sin] 35// (Barinovs et al. 
2005), but not for, e.g., [Feii]26//. Following Kaufman et al. 
(2006), it is assumed that the cross-section for H" + Si^ also ap- 
plies to fine-structure transitions of other singly ionized species. 
Concerning [Feii] (ground state ''D9/2), collisions to ^07/2 fol- 



low the above rule, but cross-sections for transitions to the next 
J are taken as 2/3, 2/4, etc. of the first one. 

Low-ionization IR lines, including dominant coolents and 
other unobserved lines, are considered in Table |6] Line identi- 
fications and observed intensities appear in Cols. 1 & 2. The re- 
sults of two computations are provided for M2 (Cols. 3-4), M2x 
(Cols. 5-6) and M2x2 (Cols. 7-8). In the first one, the excitations 
of [Sin], [Fell], [Neil], and [Am] (but not [Oi] and [Cii]) by 
collisions with H i are inhibited. In the second one, they are in- 
cluded according to the above prescription. Comparing different 
odd columns ('no') of Table |6l the rise of line intensities as the 
H I zone develops shows that the excitation of, e.g., [Si 11] by free 
electrons is still active. [Ne 11] is stable, due to the scarcity of Ne^ 
(Fig. |7}l) relative to Si^. Comparing now odd columns to even 
columns, it is seen that Hi collisions, ineffective in the 'normal' 
Hn region model M2 (except, quite interestingly, for [Neii]), 
strongly enhance the excitation rates. After subtracting emission 
from the H 11 region (M2), H i collisions contribute 75-80% of 
the excitation (virtually 100% in the case of [Ne 11]). 

Concerning [Si 11] 35//, the atomic data are not controver- 
sial and the re-calibrated flux of 68 is reasonably well defined 
(Sect. |6A2] |. Then M2x2 (Col. 8 of Table |6]l is excluded, while 
M2x (Col. 6) or an even weaker soft X-ray source (closer to the 
CHANDRA extrapolation) can account for [Si 11]. Although the re- 
markably coherent predictions for [Si 11] 35// and [Fe 11] 26// are 
partly fortuitous, they are consistent with (1) the [Fe 11] collision 
strength is correctly guessed, (2) Fe/Si is the same in the ionized 
and neutral gas, and (3) the excitation by soft X-ray heating is 
viable. Concerning [Neii] 12.8//, the discrepancy with observa- 
tion (factor 0.6) may not be significant, as the line is weak and 
its detection in the Jow-resolution mode is not taken as certain 
by Wu07. The collision strength may be too small. 

Summarizing, a plausible extrapolation to soft X-rays of 
the CHANDRA flux Can provide an explanation to the relatively 
large intensity of [Si 11] 35// and other fine-structure lines in 
IZw 18. This is considerable support to the general picture 
of photoionization as the overwhelmingly dominant cause of 
heating of the H 11 region, since heating by conversion of me- 
chanical energy appears unnecessary even in regions protected 
from ionization and heating by star radiation. Full confirmation 
should await reliable collisional excitation rates by H i for fine- 
structure lines of all singly ionized species. The soft X-rays from 
IZw 18NW have Httle effect on the Hii region: both [Oi] 6300 
and [Oiv]25.9// are increased by 3% and [Fevi] by 30%. The 
1.7% enhancement of [O iii] 4363 is of interest (Sect. | 



6.6. Stability of results 

The relative stability of the predicted line intensities is a conse- 
quence of the set of constraints (Table [T]i. Allowing for a range 
of values, a broader variety of results could be obtained. Are 
conclusions dependent on input data? 

The only basic line showing substantial variability in differ- 
ent spectroscopic studies is [0 11] 3727. This line is sometimes 
found to be stronger than the adopted value (VI98; TI05). In a 
new variant M2v of Model M2, the observed [0 11] intensity is 
assumed to be 20% larger than in Tables |3]&|4] and the covering 
factors are left unchanged. The inevitable 18% increase of the 
akeady too strong line [0 11] 7320H-30 is not too significant, con- 
sidering that the A 7325 flux is very uncertain and may not corre- 
spond to a slit position with stronger A 3727. Owing to the larger 
fractional abundance of O^, O/H is increased by 5%. The larger 
weight of low-ionization layers induces a 3% increase of Ne/H 
and a 13% decrease of N/H and Fe/H, since the [N 11] and [Fe iii] 
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Fig. 8. Cumulative mean squared relative temperature fluctua- 
tions (Peimbert, 1967) weighted by A^eXN(H+) (solid line), 
A^eXN(0+) (dotted line) and N^xH{0^^) (dashed line) versus t 
in IZw 18 NW models: (a) A^O, (b) A^l, (d) Ml; and versus R: 
(c) M2. Note the expanded vertical scale in panel (a). The short 
horizontal bars to the right of panel (d) correspond to the means 
over the full model. Vertical arrows as in Fig. [3] 



intensities were left unchanged. Both [Ariv] and [S iv] decrease 
by a few %, whilst [Oiv] increases by 6% and both [Sii] and 
[Oi] increase by ~ 11%. Ar/H and the [Sm] lines increase by 
only 1% and [Oiii]/l4363 is unchanged. 

In a more extreme example, M2cvx, with an assumed [On] 
intensity of 322 instead of 238 (factor 1.35), the /'s as in M2c 
and the SED as in M2x (reconverged = 2.5 and PiJPout - 
9), [Om] is exactly matched again, [Oi] is +2% off", [Oiv] +6%, 
[S IV] +87% and [S ii] only -8%. 

Thus, changes are moderate and turn out to alleviate diffi- 
culties noted in Sects. [63] & 16.41 e.g., the weakness of the [S ii] 
doublet. The computed r[0 iii] is robust. 

6. 7. Elemental abundances 

To first order, O/H reflects <Te{0^+)>, related to r([Oiii]), i.e., 
the predicted [Oiii]/14363 intensity. Models M3 and M2c both 
almost exactly fit /1 4363 and share the same O/H - 1.62x10"'', 
which is the best estimate, provided that (1) oxygen lines have 
been given optimal observed intensities, (2) these models faith- 
fully represent the H ii region, and (3) collision strengths are ac- 
curate. 

Concerning line intensities, [O iii] 5007 is quite stable in 
different spectra of I Zw 1 8 NW and the reasonably large, yet 
representative, ratio r([Oiii]) (Sect. 13.31 1 is taken for granted, 
since our objective is deciding whether this specific ratio can 
be consistently explained assuming photoionization by stars. In 
model M2cvx (Sect. WM . which also fits exactly the [Oiii] lines, 
[O n] 3727 was assumed to be enhanced by 35%, leading to O/H 
= 1.74x10-^ 

Concerning models, the difference between the di- 
rectly derived from r([Oiii]), re([Oiii]) = 19 850 K, and the 
NsXN{0^+) weighted average, <T^iO^+)> = 19 650 K, corre- 
sponds to a formal f^([Om]) = 0.012, similar to the computed 
j2^q2+-) _ 0.010. This difference makes only 1% difference for 
q2+^+ ^jjJ empirical estimate neglecting should nearly 
coincide with model results for this ion. A major feature of the 
Te profile is that the difference <TJP^+)> - <Te(0^)>, which 



was only 300 K in A^O and 3 200 K in A^l , is 6 600 K in best mod- 
els. Models are essential in providing a to derive O^/H^ from 
[O ii] 3727, as Tf.{[0 ii]) is poorly determined from the uncertain 
[On] 7325 and f^([Oii]) is inaccessible. 

A 'canonical' O/H for lZwl8NW is 1.46x10-^ (SK93, 
ICF99), 11% less than the present 1.62x10 ^ out of which 4% 
are due to collisional excitation of H/3 and the remaining 7% 
could be a non-trivial consequence of the relatively large f^'s 
obtained for some ions in the present models (Table |2j Fig.[8]l, 
although differences in collision strengths may also intervene at 
the 2% level (AppendixlOlli. 

The silicium and sulfur abundances were not fine-tuned 
in models. The computed Sim] flux suggests dividing the as- 
sumed Si/H by 1.25. Concerning S/H, [Sii] is underestimated, 
[S III] globally underestimated and [S iv] neatly overestimated. 
Correcting the ionization balance [S iv]/[S iii], S/H should be di- 
vided by 1.3 in best models (Sect.|6A3), but [S ii] is then under- 
estimated. Since a combination of effects may explain the over- 
estimation of [S iv], the model S/H is tentatively divided by 1.15. 
Similarly, [Ariv]/[Ariii] is best accounted for if Ar/H is divided 
by 1.13 (Sect. |63] |. but the [Ariv] line is weak and the adopted 
correction factor is 1.06. Thus, S/Ar in IZw 18 is within 10% of 
the solar value, in agreement with a conclusion of Stevenson et 
al. (1993). The iron abundance relies on the [Feiii] lines, since 
[Fe ii] 26;U does not arise from the H ii region, while the [Fe iv] 
and [Fev] intensities are uncertain. The [Feiii] intensities are 
from a spectrum in which [O ii] 3727 is stronger than on average 
(T105). In variant M2cvx (Sect. l6^ . where the intensity of [O ii] 
is multiplied by 1.35, both the predicted [Feiv]4906 and Fe/H 
are divided by 1 .4. This lower Fe/H is adopted. 

Solar abundances are tabulated by Asplund et al. (2005, 
AGS05). The compilation by Lodders (2003) is in substantial 
agreement with AGS05 (-1-0.03 dex for all O-S elements of in- 
terest here and -1-0.02 dex for Fe relative to H), except for Ar/H 
(-1-0.37 dex). The larger argon abundance is convincingly advo- 
cated by Lodders (2003). Ar/O is adopted from this reference. 
Then Ar/H coincides with the value listed by Anders & Grevesse 
(1989). The shift of O/H from Anders & Grevesse (1989) to 
AGS05 is -0.27 dex, out of which -0.07 dex corresponds to the 
change from proto-solar to solar abundances. Shifts for X/Fe are 
~ -0.20 dex for N, O, Ne, -0.1 1 dex for C, -0.07 dex for S and 
~ 0.0 for other elements of interest. 

In Table |7] the present model abundances by number 
12-Hlog(X/H) for IZw 18 NW ('M') are provided in Col. 2. The 
abundances X/O relative to oxygen from models (Col. 3, 'M') 
are compared to empirical values obtained by 1T99 (Col. 4, TT'). 
The brackets [X/Y] = log(X/Y) - log(X/Y)a from models are 
given in Cols. 5 & 6 (Y = H and Fe). Finally [X/Fe] is provided 
for Galactic Halo stars with [Fe/H] ~ -1.8 (Col. 7, 'H'). Despite 
considerable efforts to include 3D and non-LTE effects in the 
study of line formation in cool stars, astrophysical descriptions 
and atomic data may still entail uncertainties in stellar abun- 
dances {e.g., Fabbian et al. 2006), particularly for nitrogen. Also, 
a subpopulation of N-rich stars is well identified (e.g.. Carbon et 
al. 1987). 

Comparing Cols. 3 & 4, the model and empirical X/O agree 
to about 0. 1 dex. The present C/O is close to the one obtained by 
Garnett et al. (1997), who claim that C/O is anomalously large in 
IZw 18. 1T99 argue that the subregion of IZw 18 NW observed 
with the HST is especially hot according to spatially resolved 
MMT data and that C/O is therefore small. Nonetheless, 1T99 
also derive an exceedingly low O/H at the same position "be- 
cause of the higher 7^", which poses a problem of logics since 
there is a priori no link of causality between and O/H within 
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IZw 18. The [C/O] — -0.39 resulting from the present model 
is indeed marginally incompatible with the up-to-date [C/O] - 
-0.57 + 0.15 corresponding to Galactic Halo stars with [O/H] 
= -1.45 (Fabbian et al. 2006). This 'large' [C/O] is analysed 
by Garnett et al. (1997) in terms of carbon excess, suggesting 
that an old stellar population managed to produce this element, 
then challenging the view that IZw 18 is genuinely young {e.g., 
IT04), a view also challenged by Aloisi et al. (2007). From mod- 
els, [C/Fe] appears to be identical in IZw 18 and halo stars of 
similar metallicity (Cols. 6 & 7). The relatively large [C/O] in 
IZw 1 8 is due to a relatively small [O/Fe]. This is indirectly con- 
firmed by the agreement between I Zw 1 8 and halo stars for all 
elements beyond neon (argon should follow lighter or-elements). 
The [X/0]'s (= [X/H] + 1.45) are the usual basis to discuss 
elemental abundances in BCD's and nebulae. Exceptionally, in 
I Zw 18, the abundances of iron and heavy a-elements are in har- 
mony, allowing to consider the oxygen abundance with respect 
to metallicity, instead of defining metallicity by means of oxygen 
itself. Any iron locked into dust grains would further decrease 
[O/Fe] in IZw 18. Apparently, for sufliciently low metallicity of 
the ISM and/or sufficient youth of the host galaxy, iron does not 
find paths to efficiently condense into dust. Alternatively, dust 
grains may be destroyed by shocks. 

6.8. Overall evaluation of models 

Line intensities are generally well accounted for, although 
[S iv]^ 10.5;u is overpredictedby 90-100% (Sect. |6l3] ). 

The freedom left in the parameters describing the SED and 
the shell acts at the few % level upon the [Om]/l4363 predicted 
intensity. Thus, the calculated [O ni]/l4363 shifts from 96.0% of 
the observed intensity in M2 (72 = 4x10'' K, = 0.26, 0.30) 
to 99.6% in M2c = 0.22, 0.60) and 99.8% in M3 (T2 = 
5xlO''K, = 0.23, 0.50). Adding less than 1% luminosity as 
soft X-rays (e.g., M2x, compared to M2, Sect. 16. 5t , results in 
H-1.7% for [Oiii]/l4363. Also, increasing He/H from 0.080 to 
the possibly more realistic value 0.084 (Peimbert et al. 2007), 
[Om]/i4363 is enhanced by a further -1-0.6%. Since both the X- 
ray and He/H corrections are more than plausible, it is relatively 
easy to reach 100-102% of the observed [Ora]/l4363 intensity 
in the assumed configurations. 

However, models tell us that r([Om]) can hardly be larger 
than the observed value. It may prove necessary to consider alter- 
native gas distributions, e.g., in case the [S iv] misfit is confirmed 
by more accurate observational and theoretical data. Assuming 
larger densities in the difiiise medium (Sect. 16.4.3b and/or con- 
sidering thick filaments closer to the source tend to penalize 
[Om]^4363. 

At this point, uncertainties on [Oiii] collision strengths Q 
need be considered (Appendix lD.2l l. Using fi's by Lennon & 
Burke (1994, LB94) instead of Aggarwal (1993, Ag93), the 
computed /1 4363 would be 2.1% smaller and more difficult to 
explain. On the other hand, using Q(^P- 'D) from Ag93 and 
Q(^P - 'S) from LB94 would enhance all computed /1 4363 in- 
tensities by 2.1%. Concerning transition ''P - 'S which controls 
.14363, both Ag93 and LB94 find a 5-6% increase of Q from 
2x10''K to 3x10^ K, suggesting the influence of resonances. If 
for some reason the energy of these resonances could be shifted 
down, there would be room for a few % increase of Q at 2x10"* K 
compared to the current value, then introducing more flexibility 
in the present model of IZw 18. 

Summarizing, the hypothesis of pure photoionization by 
stars in the form explored here is perfectly tractable, but the 



models approach a limit. This is in a sense satisfactory, consid- 
ering that IZw 18 is an extreme object among BCD's, but suffi- 
cient flexibility in choosing solutions is worthwile. An analysis 
of how the computed r([Oiii]) can be influenced shows that, in 
the case of IZw 18, possible variations of 'astrophysical' ori- 
gin are of the same order as the uncertainties affecting the Q's. 
Since the set of computed r([Oiii]) tends to be down by 2-3% 
relative to observation, it is legitimate to question the Q's. Now, 
the recent re-evaluation of the distance to IZw 18 by Aloisi et 
al. (2007) may offer an 'astrophysical alternative': multiplying 
D, Ri and Rf by 2'^^ and the luminosty by 2, the relative volume 
increase leads to smaller Pi,, and t^, and, after reconvergence, 
[Oiii]/14363 is enhanced by -i-2.4%. Noneteless, [Siv] 10.5ju is 
enhanced too. 

7. Concluding remarks 

Owing to its small heavy element content, I Zw 1 8 stands at the 
high-Te boundary of photoionized nebulae. Where ionization 
and temperature are sufficiently high, the cooling is little de- 
pendent on conditions, except through the concentration of H", 
controlled by density. Therefore, in the photoionization model 
logics, Tf. is then a density indicator, in the same way it is an 
O/H indicator in usual Hii regions. It is for not having recog- 
nized implications of this new logics, that low-metallicity BCD 
models failed. 

In a photoionization model study of IZw 18NW, SS99 em- 
ployed a filling-factor description and concluded that re([Oiii]) 
was fundamentally unaccountable. The vogue for this simple 
description of the ionized gas distribution resides in its appar- 
ent success for usual Hn regions, a success falling in fact to 
the strong dependence of cooling on abundances. Universally 
adopted along past decades, this concept led all authors to con- 
clude that photoionization by hot stars did not provide enough 
energy to low-Z GEHIIRs. This conclusion is in line with a 
movement of calling into question photoionization by stars as 
the overwhelmingly dominant source of heat and ionization in 
gaseous nebulae, a movement cristalizing on the 'f^ problem' 
(Esteban et al. 2002; Peimbert et al. 2004), since the presence of 
Te fluctuations supposedly Jarger than those reachable assuming 
photoionization by stars implies additional heating. 

A conclusion of the present study is that the gas distribu- 
tion is no less critical than the radiation source in determin- 
ing the line spectrum of Hn regions. Assuming pure photoion- 
ization by stars, the remarkable piece of information carried 
by the large re([Oiii]) of IZwl8NW is that the mean den- 
sity of the [Om] emitting region is much less than A^e([Sii]), 
a low A^e confirmed by line ratios [Oiv]/125.9yL//Heii/l4686 
and [Fem]^4986/[Feiii]^4658. IZwl8NW models compris- 
ing a plausible SED and respecting geometrical constraints can 
closely match almost all observed lines from UV to IR, includ- 
ing the crucial [O m]/} 4363 ([S iv]A IQ.Sfi is a factor 2 off, how- 
ever). Thus, extra heating by, e.g., dissipation of mechanical en- 
ergy in the photoionized gas of low-metallicity BCD galaxies 
like IZw 18 is not required to solve the 're([Oiii]) problem'. 
Moreover, since low-ionization fine-structure lines can be ex- 
plained by soft X-rays, (hydrodynamical) heating is not even re- 
quired in warm H i regions protected from ionization and heating 
by star radiation. 

As a final note, on close scrutiny, the solutions found here 
tend to be just marginally consistent with observed r([Om]). 
Given the claimed accuracy in the different fields of physics and 
astrophysics involved, postulating a mechanical source of heat- 
ing is premature, whereas a 2-3% upward correction to the col- 
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lision strength for transition 0^+(^P- 'S) at ~ 2x10"^ K is an 
alternative worth exploring by atomic physics. Yet, another pos- 
sibility is a substantial increase of the distance to 1 Zw 18. Owing 
to accurate spectroscopy and peculiar conditions in IZw 18, im- 
portant astrophysical developments are at stakes in the 5% un- 
certainty attached to [Om] collision strengths. 

If photoionized nebulae are shaped by shocks and other hy- 
drodynamic effects, this does not imply that the emission-line 
intensities are detectably influenced by the thermal energy de- 
posited by these processes. Unravelling this extra thermal en- 
ergy by means of spectroscopic diagnostics and models is an 
exciting prospect, whose success depends on a recognition of all 
resources of the photoionization paradigm. Adopting the view 
that photoionization by radiation from young hot stars, includ- 
ing WR stars, is the only excitation source of nebular spectra in 
BCD galaxies, yet without undue simplifications, may well be 
a way to help progresses in the mysteries of stellar evolution, 
stellar atmosphere structure, stellar supercluster properties, gi- 
ant H n region structure and, at last, possible extra sources of 
thermal energy in BCDs. 
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Appendix A: Models for other GEHIIRs 

Models of GEHIIRs include studies of individual objects and 
evolutionary sequences for large samples. Examples ordered by 
decreasing O/H are first reviewed. 

A.1. Individual GEHIIRs 
A.1.1. O/H > 1.5x10 ■* 

Gonzalez Delgado & Perez (2000) successfully model NGC 604 
(O/H = 3x10""*) in M33 as a radiation-bounded sphere (ra- 
dius 20-llOpc) of density 30 cm"^ and filling e = 0.1, both 
[Oi]/I6300 and [Oiii]/l 4363 being explained. 

Garcia- Vargas et al. (1997) successfully model cir- 
cumnuclear GEHIIRs in NGC 7714 as thin, constant- 
density (A^H ~ 200 cm"''), radiation-bounded shells with 
O/H = (2-3)xl0"^: r([Oiii]) is accounted for within errors and 
[O i] is just moderately underestimated. The nuclear GEHIIR of 
NGC 7714 is modelled by Gonzalez Delgado et al. (1999) as a 
full sphere with very small e. The adopted O/H = 3x10 is too 
large since [Oiii]/l4363 is underpredicted. Obviously, a better 
fit to the available optical line spectrum could again be achieved 
for the nucleus. Gonzalez Delgado et al. (1999) are probably not 
well founded to invoke extra heating by shocks. 

Luridiana & Peimbert (2001, LPOl) propose a photoion- 
ization model for NGC 5461 (O/H = Z5xlO **), a GEHIIR in 
MlOl. As for NGC 2363 (AppendixjAlllll, LPOl apply an aper- 
ture correction to their spherical model. The H/3 and r([S ii]) spa- 
tial profiles are reproduced with a Gaussian density distribution 
of very small e and high inner density - 500 cm"^ compared to 
A^e([S ii]) ~ 150 cm"^ - meant to increase the inner fraction. 



In this way, the /1 5007, /1 4363 and /1 3727 fluxes restricted to 
the theoretical slit can be accounted fofl but not [Oi], "a not 
unusual fact", nor [S ii], which, unlike a belief of LPOl, is not 
enhanced by increasing the primary flux below 1 .0 ryd. As noted 
by LPOl, the outputs of their model are strongly dependent on 
the density structure. 

From their sophisticated study of NGC 588 (O/H = 2x10"'*), 
Jamet et al. (2005, JS05) conclude that "the energy bal- 
ance remains unexplained". This negative conclusion is es- 
sentially based on the fact that re([Om]), from the ratio 
/1 4363 A/^ 5007 A, is observed to be larger than re([Oiii], IR), 
from ^5007 A//188yum, by ATob, = 2 700+700 K, while the cor- 
responding AT, nod is only 1 400+200 K in models, which are 
otherwise satisfactory, accounting reasonably well for re([Oiii]) 
and the distribution of ionization (models DDI and DDH ex- 
hibited by JS05, who carefully consider uncertainties related to 
the SED and the small-scale gas distribution). Considering the 
difficulty of calibrating the ISO-LWS fluxes relative to the opti- 
cal and the sensitivity of the [Oiii]/l88//m emissivity to A^e, the 
400 K gap between hJobs and hTmod is not a sound basis to claim 
the existence of an energy problem. If diagnostics based on IR 
lines are desirable, the energy problem raised so far in GEHIIR 
studies is not related to these lines. Instead of a heating problem 
as in IZw 18, the model presented by JS05 could be facing a 
cooJing problem, since the computed re([Oiii], IR) is too high. 

An additional energy source is not clearly needed in the cases 
of NGC 588, NGC 5461, NGC 7714 and NGC 604. 

A. 1.2. O/H < 1.5x10-4 

Relafio et al. (2002) provide an inventory of NGC 346, a GEHIIR 
of the SMC (O/H - 1.3x10 **). Their spherical, constant-density, 
matter-bounded photoionization model, whose only free param- 
eter is a filling factor e (alia SS99), accounts for the escape of 
ionizing photons, but underpredicts collisional lines, especially 
[O III]/! 4363. After unsuccessful variations on geometry, the au- 
thors preconize, following SS99, an additional source of energy. 

After an extensive exploration of photoionization models 
with filling factor for the bright GEHIIR NGC 2363 (O/H = 
8x10"^), Luridiana et al. (1999, LPL99) conclude that they can- 
not find a solution unless they introduce fluctuations by hand, 
i.e., they assume a larger than the one intrinsic to their model. 
This f^, intended to enhance the computed [Oiii]/l4363, is jus- 
tified by the fact that the observed Paschen jump temperature 
is less than re([Oiii]) and supported by a self-consistency argu- 
ment: a larger leads to a larger O/H, hence a larger number 
of WR stars, hence (1) a larger injection of mechanical energy, 
supposed to feed the temperature fluctuations themselves, and 
(2) a larger photon flux above the He^ ionization limit, useful to 
increase Heii/14686. However, as acknowledged by Luridiana 
et al. (2001), the WR star winds make a poor job to generate 
a significant in NGC 2363. Also, present views suggest that 
arguments based on WR stars in low-Z galaxies were illusory 
just a few years ago (e.g., Leitherer 2006; AppendixICll. Finally, 
r([Oiii]), underestimated by only ~ 12% in the 'standard' low-Z 
model by LPL99, is divided by 2 on using the larger O/H, so that 
a relatively minor difficulty is first made much worse and then 



^ LPOl state a priori that [Om]/i4363 "almost surely has a contribu- 
tion from processes other than photoionization" and, consistently, con- 
clude that their model "fails to reproduce the observed [O m]A 4363 in- 
tensity", but both statements seem to be refuted by pieces of evidence 
they present. 
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solved by means of a f^. The slit correction advocated by LPL99 
is considered in Appendix |B.2| 

Luridiana et al. (2003, LPPC03) consider a spherical model 
for a GEHIIR of SBS 0335-052 (O/H = 2x10-^). A Gaussian 
distribution with large maximum density and small e proves un- 
satisfactory. LPPC03 then consider a 10-shell model (over 50 
free parameters, most of which pre-defined), in which each shell 
is radiation bounded and is characterized by a covering factor 
Although each shell is still given an e, the new model is equiva- 
lent to a collection of geometrically thin radiation-bounded sec- 
tors at different distances from the source (see also Giammanco 
et al. 2004) and "gracefully reproduces the constancy of the ion- 
ization degree along the diameter of the nebula". Hence, the au- 
thors are forced by observational evidences to implicitly aban- 
don the classical filling-factor approach. Nonetheless, whatever 
the complexity of these models, all of them fail to account for 
the high re([Oiii]). 

LPPC03 consider a Gaussian model for IZw 18 SE (O/H = 
1.7x10"^) with again a relatively large maximum density and, 
unlike for SBS 0335-052, a relatively large e, resulting in a rather 
compact model nebula, in which the computed re([Oiii]) com- 
pares quite well with the observed one. Unlike for the NW, the 
HST image (Cannon et al. 2002) of the younger SE Hii region 
does not show a shell surrounding an MSG. Nonetheless, consid- 
ering the strong output of mechanical energy from massive stars, 
it is likely that inner cavities already developped. The strong in- 
direct evidence for too compact a gas distribution in the model 
by LPPC03 is the notable weakness of the computed intensity of 
[O ii] and other low-ionization lines. Adopting a more expanded 
structure in order to increase [O ii], yet keeping the general trend 
of the gas distribution, the computed re([Oiii]) would be forsee- 
ably lower than in the model by LPPC03. 

A.2. Individual GEHIIRs: discussion 

LPPC03 describe the '7'e([Oiii]) problem' they face in their 
study of SBS 0335-052 ( Appendix lA. 1 .2b as "a systematic fea- 
ture" of Hii region models and, following SS99, they state 
that this problem "can be ascribed to an additional energy 
source acting in photoionization regions, other than photoion- 
ization itself". Nevertheless, re([Oiii]) seems to be account- 
able in existing photoionization models for GEHIIRs with, say, 
O/H > 1.5x10""* (Appendix lA. 1.1b . Similarly, the computed 
[Om]/14363 is correct, possibly even too large, for Hii regions 
of the LMC (Oey et al. 2000). If, despite apparent complemen- 
tarities, the re([Oiii]) and problems have different origins 
(Sect. [Til, no 'systematic feature' can be invoked. 

In modelling objects with near solar abundances, 
[Oiii]/14363 is controlled by O/H, [Oiii]/15007 by the 
'color temperature' of the ionizing radiation, [Oii]/l3727 by 
the ionization parameter, while [O i] is maximized in radiation- 
bounded conditions. For these objects, assuming a 'large' 
(constant) density, e.g., ~ A^e([Sii]), associated to an ad hoc 
e <sc 1, is often successful, although this does not prejudge of the 
relevance of the model found. Indeed, this assumption proves 
to be at the heart of the re([Oiii]) problem met in low-Z BCDs 
(Appendix IB.lb . 

A.3. Extensive analyses of BCDs 

The conclusion of an early extensive analysis based on radiation- 
bounded, low-density full sphere models for low-Z BCDs 
(Stasiriska & Leitherer 1996) is optimistic concerning re([0 iii]). 



whereas [Oi] is then qualitatively explained in terms of shock 
heating. Nevertheless, in an extension of this study to large- 
Z objects with no measured re([Oiii]), Stasinska et al. (2001) 
reinforce the energy problem raised by SS99 (Sect. 12. 2b when 
they conclude that "a purely 'stellar' solution seems now clearly 
excluded for the problem of [O iii]/H/? versus [O ii]/Hy6 as well 
as [S ii]/H/3", while, conversely, they still endorse the unproved 
statement of SS99 (Sect. 1231 ) that "strong [O i] emission is easily 
produced by photoionization models in dense filaments". 

The sequence of photoionization models proposed by 
Stasiriska & Izotov (2003, SI03) for a large sample of low-Z 
BCDs (divided in three abundance bins) illustrates views ex- 
pressed after the failure of models acknowledged by SS99 for 
IZw 18 (Sect. |2]i. In the description by SI03, an evolving syn- 
thetic stellar cluster (10^ Mq, instantaneous burst) photoionizes 
a spherical shell of constant density A^h = 10^ cm"^ at the bound- 
ary of an adiabatically expanding hot bubble. With suitable bub- 
ble properties, underlying old stellar population, aperture cor- 
rection and time evolution of the covering factor, the range of 
H/3 equivalent width (EW(H/?)) and the trends of [Om] 5007, 
[O II] 3727, [O i] 6300 versus EW(H/3) can be reproduced for the 
high-Z bin (O/H ~ 1.5x10 "*) within the scatter of the data. 

Applying similar prescriptions to the intermediate-Z bin, the 
oxygen lines and He iiA 4686 (He ii was just fair in the first bin) 
are underpredicted. SI03 diagnose an insufficient average energy 
per absorbed photon and assume that the stellar cluster is supple- 
mented by a strong Kfi K bremsstrahlung-like radiation source, 
which solves the He ii problem (He^ is further ionized by extra 
4-5 ryd photons; see, however, AppendixICb and alleviates the 
[O i] problem (the soft X-rays further heat and widen the ioniza- 
tion front), but barely improves [Oiii] and [On]. Agreement of 
the model sequence with observation is finally restored by sup- 
posing in addition that the shell includes a time-variable oxygen- 
rich gas component attributed to self-enrichment: in the example 
shown by SI03, this component is 4-fold enriched in CNO, etc. 
relative to the original abundance and encompasses half of the 
shell mass after a few Myrs, so that one generation of stars pro- 
duced a 2.5-fold enhancement of the average abundance in the 
photoionized gas. 

This description essentially applies to the low metallicity 
bin (O/H ~ 2x10"^) of particular concern for IZw 18, but with 
even more extreme properties for the O-rich component, since it 
should be overabundant by 1 dex, resulting in a 5-fold enhance- 
ment of the final average abundance. 

A.4. Extensive analyses of BCDs: discussion 

The time scale of 0.5 Myrs for the growth of the O-rich com- 
ponent in the description by SI03 cannot directly fit in the self- 
pollution scenario since it is shorter than the stellar evolution 
time scale. Also, a sudden oxygen self -pollution of the gas is not 
observed in supernova remnants. 

The assumed X-ray power is ~ 10% of the cluster luminos- 
ity or ~2dex times the estimated X-ray rosat power (0.07-2.4 
keV) of the hot bubble fed by stellar winds and supernovae 
around a usual MSG (Strickland & Stevens 1999; Cervino et 
al. 2002). Moreover, the hot gas is generally raised at severaJ 
10^ K (Stevens & Strickland 1998). Adopting a larger tempera- 
ture, the X-ray power should be even larger, as only the softer 
radiation interacts usefully with the ionized gas. IZw 18 itself 
is a rather strong X-ray emitter in the 0.5-10 keV range, yet 
20 times weaker than the source assumed by SI03 (Thuan et al. 
2004; Sect.|63]l. 
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Apart from these problems, SI03 do not address the question 
of the intensity of [Om]/i4363. The narrow radiation-bounded 
shell adopted by SI03 usefully favours [Oi]/i6300, but makes 
the computed intensity of [O iii]/l4363 even worse than the one 
obtained by, e.g., SS99 (Sect. |2|. Moreover, adding the promi- 
nent O-rich component advocated by SI03 will (1) decrease 
[Oiii]/l4363 by a further 30-40% on average and (2) conflict 
with the existence of very low-Z BCDs, since any of them will 
be condamned to shift to the intermediate class defined by SI03 
after just 1 or2Myrs. 

If what SI03 qualify as "appealing explanations" presents 
any character of necessity for BCD models, then the hypothe- 
sis of photoionization by stars, which was already given a rough 
handling by SS99 in their analysis of I Zw 1 8, should be consid- 
ered as definitively hurried for the whole class of low-Z BCDs. 
The fact that SI03 discard [Oiii]/l4363 in their analysis con- 
firms that they endorse and reinforce views expressed by SS99 or 
Stasinska et al. (2001) and give up explaining re([Oiii]) in low- 
metallicity GEHIIRs by means of stellar radiation. However, the 
same restrictive assumption as for individual GEHllR studies 
(Appendices lA. H4A.2I ) bears on the gas distribution adopted by 
S103, since their (geometrically thin) 'high' constant-A^n model 
sphere is nothing but a zero-order approximation to a model shell 
with classical filling factor (Appendix IB.lb . 

Appendix B: On the gas distribution in GEHIIRs 

B.1. Misadventures of the filling factor concept 

For the sake of reproducing the Ha surface brightness of Hii 
regions, asssuming a gas density much larger than <Ne^>^^^, the 
'filling factor paradigm' posits that the emitting gas belongs to 
optically thin, 'infinitesimal' clumps, filling altogether a fraction 
e 1 of the volume. 

Given that the stellar evolution timescale exceeds the sound 
crossing time of H ii regions, small optically thin ionized clumps 
will have time to expand and merge into finite-size structures. If 
these structures have the original density, they are likely to have 
finite or large optical depths, in contradiction with the filling fac- 
tor concept. 

Filamentary structures, ubiquitous in Ha images of nearby 
GEHIIRs, are often taken as justifications for introducing e in 
photoionization models. However, (1) the geometrical thickness 
of observed filaments is consistent with radiation-bounded struc- 
tures and (2) an individual filament most often emits both high 
and low ionization lines {e.g., Tsamis & Pequignot 2005). The 
filling-factor description is flawed. 

A GEHIIR may well be a collection of radiation-bounded 
filaments embedded in coronal and photoionized diffuse media. 
The idea behind assuming this configuration is that only the ion- 
ized 'atmospheres' of long-lived, radiation-bounded, evaporat- 
ing structures will possibly maintain a substantial overpressure 
relative to their surroundings (a similar idea applies to the "pro- 
plyds" found in Orion; e.g., Henney & O'Dell 1999). 

The filling factor concept fails on both theoretical and ob- 
servational grounds. Nevertheless, introduced as a technical tool 
to manage A^e diagnostics like r([S ii]), e came to be improperly 
used to adjust the local ionization equilibrium of the gas through 
A^H^ in ™ effort to overcome problems of ion stratification gen- 
erated by the filling factor description itself (Appendix lB.2l i. 

The re([Oiii]) problem met in oxygen-poor GEHIIRs may 
relate to the loss of plasticity affecting photoionization mod- 
els, as the dependence of gas cooling on abundances van- 
ishes. Then, cooling depends on the relative concentration of 



H" (collisional excitation of Lya), controlled by the local Nu- 
Hence, the (improper) freedom on e is eroded. Moreover, if 
density is not uniform, A^e([Sn]) is a biased estimate for A^e 
in the bulk of the emitting gas, since ions will belong to 
dense, optically thicker clumps. Emission from an interclump 
medium with A^e <A^e([Sii]) will selectively enhance the com- 
puted [O III]/! 4363 intensity. LPL99 state promisingly that their 
model includes "denser condensations uniformly distributed in 
a more tenuous gas", but in practice only the condensations 
emit. This restriction is shared by virtually all published mod- 
els for low-Z GEHIIRs. While the assumed density of the emit- 
ting gas can be orders of magnitude larger than <Ne^>^^^, emis- 
sion from a lower density gas is neglected by construction (The 
study by JS05 is an exception, but NGC 588 is not low-Z; 
Appendix lA.l.ll i. The re([0 iii]) problem suggests lifting this re- 
striction. 

B.2. Spheres, slits, filling factor and stratification 

Spherical models raise the question of how to compare computed 
spectra with nebular spectra observed through, e.g., a narrow 
slit. LPL99 advocate extracting emission from that part of the 
sphere which would project on the slit. Despite obvious prob- 
lems with non-sphericity, LPL99 and others argue that this pro- 
cedure would at least allow weighting the contributions from 
low- and high-ionization zones in a more realistic manner Using 
the classical filling factor concept (Appendix IB.lb in GEHIIR 
models, ion stratification spreads over the whole nebula and the 
[Oi] emission is eff'ectively confined to outer layers, in which 
the primary radiation eventually vanishes. If, on the contrary, the 
emitting gas belongs to radiation-bounded filaments distributed 
within the nebula, then ion stratification disappears to first or- 
der. Radial ionization gradients, if any, are no more related to 
a progressive destruction of primary photons along the full ra- 
dial extension of the nebula, but to changes in (local) average 
ionization parameter. 

LPL99 conclude that [Oi] is due to shock excitation in 
NGC 2363 (Sect. lATL2b because the computed intensity is weak 
in their theoretical slit extraction. Nonetheless, the [O i] intensity 
is fairly correct in their global spectrum. This apparent failure of 
their photoionization model may well be due to the unfortunate 
combination of (1) a very small e and (2) the extraction of a slit 
shorter than the diameter of the model sphere. Along the same 
line, LPPC03 are confronted to undesirable consequences of the 
filling factor assumption on the variation of ionization along a 
sHt crossing SBS 0335-052 (Appendix lAT2] l. 

If a geometrically defined model can hardly provide an ap- 
proximation to a complex H ii region, thus casting doubts on the- 
oretical slit extractions, global spectra are less sensitive to geom- 
etry, owing to conservation laws. 

Moreover, in computing 1-D photoionization models, the 
(spherical) symmetry enters only in the treatment of the diffuse 
ionizing radiation field, which is generally not dominant in the 
total field. The diffuse field, most effective just above the ioniza- 
tion limits of H, He and He^, is relatively local at these photon 
energies (in accordance with the 'Case B' approximation) and 
little dependent on global geometry. Let us define an "elemen- 
tary spherical model" (for given SED) by a radial density distri- 
bution of whatever complexity. Since the local state of the gas 
is chiefly related to the primary (radial) radiation, a composite 
model made of a judicious combination of elementary spheri- 
cal models, each of them restricted to a sector characterized by a 
covering factor, can provide topologically significant and numer- 
ically accurate descriptions of global spectra for nebulae with 
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complex structures. Defining a 'topology' as a particular set of 
spherical models with their attached covering factors, any given 
topology is in one-to-one correspondence with a global spectrum 
and a full class of geometries, since any sector can be replaced 
by an arbitrary set of subsectors, provided that the sum of the 
covering factors of these subsectors is conserved. 

Thus, a good modelling strategy for a GEHllR is one in 
which a global (probably composite) model spectrum is com- 
pared to the observed global spectrum. If only one slit obser- 
vation is available, given that the ion stratification tends to be 
relatively loose and erratic in GEHllRs, it is wise to directly use 
this spectrum as the average spectrum (together with scaling by 
the absolute H/3 flux), with the understanding that the resulting 
photoionization model will represent a 'weighted average' of the 
real object. For many practical purposes, this weighting may not 
significantly impact on the inferences made from the model, un- 
less the slit position is exceedingly unrepresentative. 



Appendix C: He lU 4686, WR stars and SEDs 

IZw 18 harbours Wolf-Rayet (WR) stars (Legrand et al. 1997; 
Izotov et al. 1997; de Mello et al. 1998; Brown et al. 2002). 
WR stars have been challenged as the sole/main cause of nebu- 
lar He iiA 4686 in BCDs on the basis of a lack of correlation be- 
tween the occurence of this line and the broad 'WR bumps' {e.g., 
Guseva et al. 2000). The study of WR stars is experiencing a rev- 
olution (Maeder et al. 2005; Meynet & Maeder 2005; Grafener 
& Hamann 2005; Vink & de Koter 2005; Crowther 2007) after 
the realization that (1) rotation of massive stars favours enhanced 
equatorial mass loss, element mixing by shears, and angular mo- 
mentum transport by meridian circulation, (2) low-Z massive 
stars tend to be fast rotators and accelerate as they evolve off the 
main sequence, so that the lower mass limit for a star to become 
a WNE star is much reduced, and (3) for a given type of WR star, 
the mass loss is lower for lower metallicity (Fe/H, not O/H), with 
three consequences: the broad WR features are less evident for 
low metallicity (weaker optical continuum and smaller EW of 
WR bumps), the duration of the WR stage can be longer, and 
the EUV luminosity is larger due to reduced blanketing effect. 
Thus, the above lack of correlation can now be partly ascribed 
to a bias, related to the tendency of WR star atmospheres to dis- 
play less prominent optical signature when they emit more EUV 
radiation. The WR star population of IZw 18 and the ability of 
these stars to emit radiation beyond 4ryd have almost certainly 
been grossly underestimated (Crowther & Hadfield 2006). 

Other observations, e.g., for SBS 0335 052E (Izotov et al. 
2001b; Izotov et al. 2006b) are still taken as evidence for Hen 
excitation by radiation from very fast shocks: (1) the He ii line is 
broader than other nebular lines, (2) the He ii emission is spread 
out far away from the main MSCs, and (3) is larger in He ii 
emitting area, hence at large distances from the main ionizing 
sources. These findings are definitely no compelling arguments 
against photoionization by WR stars. The larger He ii line width 
indicates larger turbulence and/or velocity gradients, not neces- 
sarily shocks. That Tg is observed to be larger in Hen emitting 
gas is in agreement with photoionization models. The spatial ex- 
tent of He II may reflect the distribution of a few WR stars, which 
may not belong to the main cluster and may not be easily de- 
tected (Crowther & Hadfield 2006). Alternatively, He ii can be 
produced far from the ionizing stars if the medium is porous and 
permeated by low density, optically thin gas, e.g., along a galac- 
tic wind outflow (Izotov et al. 2006b). The picture of a galactic 
wind also suggests an explanation for the He ii width. 



Photoionization models are test beds for ionizing radiation 
sources, but inferences on the physics of GEHllRs should not 
depend on uncertain SEDs. Existing synthetic star clusters are 
inadequate to model IZw 18. Apart from known problems with 
star sampling (Cervino et al. 2003; Cervino & Luridiana 2006), 
limited knowledge of the history of actual MSCs and current un- 
certainties about WR stars, new free parameters (initial angular 
momentum and magnetic field of individual stars; rate of bina- 
rity) will broaden the range of possible SED evolutions, while 
collective effects in a compact cluster of massive stars may in- 
fluence the output of ionizing radiation far from it, due to high- 
density stellar winds (Thompson et al. 2006). 

These comments justify (1) the assumption of an excitation 
of He II solely by WR stars and (2) the use of a flexible analytical 
SED for I Zw 18 NW (Sect.gjJ. 

Appendix D: Atomic data 

D.1. Collisional excitation of H i 

Collision strengths Q.{ls-nl) {n < 6;l < n) for Hi are taken 
from Anderson et al. (2000, ABBSOO). The O's for ls-2s and 
1 s-2p are much larger than for the next transitions 1 s-nl and are 
not controversial. The main coolent agent in low-Z BCDs should 
be correctly implemented in all codes. Nonetheless, in the con- 
ditions of IZw 18, the results for transitions 1-2 by ABBSOO 
are about 10% larger than those carefully fitted by Callaway 
(1994), giving an estimate of possible uncertainties. The adopted 
data tend to enhance the cooling with respect to earlier data and 
to (conservatively) worsen the 're([Oiii]) problem'. Total Q(l- 
«)'s listed by Przybilla & Butler (2004) virtually coincide with 
ABBSOO values for 1-2, confirming the Hi cooling rate, but di- 
verge from ABBSOO for n > 2 and increasing T^. similarly to 
early, probably wrong, data (see Pequignot & Tsamis, 2005). 

D.2. Collisional excitation of[0 m] 

Effective collision strengths O obtained over past 50 years are 
listed in Table IDT] at four Te's for transitions ^P- 'D, "S 
and 'D - 'S. Aggarwal & Keenan (1999) did not feel it neces- 
sary to update earlier values by Aggarwal (1993; Ag93), almost 
contemporary with Lennon & Burke (1994). The ratios of the re- 
cent values are given in Table ID. 11 The differences are over 4% 
for 3p - 'D and 10% for 'D - 'S (6% in 1 Zw 18 conditions), but 
the latter has no influence at low A^e- nebu includes a fit better 
than 0.5% to Ag93 data. 

The [Oiii] transition probabilities used in nebu are from 
Galavis et al. (1997, GMZ97). The accuracy of the Opacity 
Project (OP) data for these transitions is 8-10% (Wiese et al. 
1996). Coherently, the much more elaborate results by GMZ97 
diff^er from the OP results by 9.6% and 5.5% for A(iD - 'S) and 
A('P - 'S) respectively. Would A('D - 'S) change by as much as 
5%, the branching ratio of [O iii]A 4363 would change by 0.6%. 

Thus, discrepancies not exceeding 5% exist among different 
calculations (3% for Q ratios), suggesting that uncertainties on 
the computed r([Oiii]) are probably < 5%. The 25-30% under- 
estimation found by SS99 is not due to erroneous atomic data. 

D.3. Miscellaneous data 

The adopted table for radiative and dielectronic recombina- 
tions is limited to the 11 sequences H-like-Na-like (Badnell 
2006). Dielectronic rates for [S ii]-[S iv] are obtained by Badnell 
(1991), but total recombination coefficients for (recombined 
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ions) [Sin], [Sii], [Sm], [Arv], [Fen]-[Fev], are taken from 
Nahar and co-workers (Nahar, 2000 and references cited). The 
recombination rate for [S n] used in this and previous nebu com- 
putations is 1.15 times the Nahar's value. Empirical total rate 
coefficients based on PN models (Pequignot, unpublished), im- 
plemented in NEBU for a decade, are 5 and 8 times the radiative 
ones for [Ar ii] and [Ar iii] respectively. A larger factor is sus- 
pected for [Arm] at high Tg. 

Collision strengths of special mention include those for 
[On] (Pradhan et al. 2006; also Tayal 2006b), [Oiv] (Tayal 
2006), [Siii] (Tayal & Gupta 1999), [Siv] (Tayal 2000) and 
[Fev] (Woste et al. 2002). Collisions with H*' are considered 
in Sect. 16.51 Charge exchange rates with H'^ for O^^ and N^^ are 
now from Barragan et al. (2006). 



Table 3. Model outputs for I Zw 18 NW: A^O and Ml 



Line id./Models 


Obs. 


NO 


NO/O 


M2 


M2/0 


Absolute fluxes (1(H8) in 


10 ''* erg cm - s ' 


) 




Km 


5.6 


5.6 


1.00 


5.6 


1.00 


1.43 GHz/mJy 


0.433 


0.240 


0.55 


0.256 


0.59 


8.45 GHz /mJy 


0.286 


0.202 


0.71 


0.216 


0.76 


Relative line fluxes (wavelengths in A or iim) 




Hi 4861 


1000. 


1000. 


1.00 


1000. 


1.00 


Hi 6563 


2860. 


2840. 


0.99 


2830. 


0.99 


Hi 4340 


461. 


473. 


1.03 


473. 


1.03 


Hi 4102 


266. 


267. 


1.00 


267. 


1.00 


Hi 1215 (/lO) 




2950. 




3010. 


- 


ni2hv(l\Q) 




1550. 




1620. 


- 


He: 3888 


90.4 


89.6 


0.99 


90.5 


1.00 


He 1 4471 


21.4 


34.9 


1.63 


35.2 


1.64 


He: 5876 


61.1 


92.0 


1.36 


91.3 


1.35 


He: 6678 


25.3 


25.3 


1.00 


25.6 


1.01 


He: 7065 


24.4 


23.4 


0.96 


22.9 


0.94 


He: 10830 




251. 




190. 


- 


Hen 4686 


36.8 


36.8 


1.00 


36.8 


1.00 


Cm] 1909+07 


467. 


467. 


1.00 


467. 


1.00 


Sim] 1882+92 


270.: 


229. 


0.85 


340. 


1.26 


Aim 1855+63 


111.: 


42.9 


0.39 


82.9 


0.75 


Oiii] 1664 


<230 


127. 


>.5 


208. 


>.9 


Civ 1549 


510.: 


74.2 


0.14 


334. 


0.65 


Si IV 1397 )° 


<300 


22.7 


>.l 


127. 


>.5 


Oiv] 1398 ) 




2.0 




24.9 


* 


[Nil] 6584+48 


9.2 


9.2 


1.00 


9.2 


1.00 


[Oi] 6300+63 


8.5 


0.12 


0.01 


8.6 


1.01 


[On] 3726+29 


238. 


238. 


1.00 


238. 


1.00 


[On] 7320+30 


6.3: 


9.5 


1.50 


7.5 


1.18 


[Om] 5007+.. 


2683. 


2680. 


1.00 


2680. 


1.00 


[Oiii] 4363 


65.9 


47.8 


0.73 


63.2 


0.96 


[Om] 51.8 yum 




174. 




137. 


- 


[Om] 88.3 yum 




216. 




213. 


- 


[Oiv] 25.9 /im 


49.1 


18.2 


0.37 


47.8 


0.97 


[Neil] 12.8 pm 


9.0: 


1.3 


0.14 


1.9 


0.21 


[Nem] 3868+.. 


191. 


191. 


1.00 


191. 


1.00 


[Neiii] 15.5 yum 


45.7 


60.0 


1.31 


48.9 


1.07 


Mgi] 4571+62 


<3.0 


1.5 


>.5 


1.2 


>.4 


[Sin] 34.8 yum 


157. 


4.7 


0.03 


22.0 


0.14 


[Sii] 6716 


22.5 


6.7 


0.30 


17.6 


0.78 


[Sii] 6731 


16.9 


5.1 


0.30 


13.1 


0.78 


[Sii] 4068 ) 


3.7 


1.1 


0.41 


2.2 


0.99 


[Fev] 4071 ) 


* 


0.4 




1.5 


* 


[Sm] 9531+.. 


114. 


130. 


1.15 


113. 


0.99 


[Sm] 6312 


6.7 


6.0 


0.89 


5.7 


0.85 


[Sm] 18.7 yum 


28.0 


32.2 


1.15 


26.2 


0.95 


[Sm] 33.5 yum 


120. 


54.5 


0.45 


48.0 


0.40 


[Siv] 10.5 yum 


48.0 


41.7 


0.87 


92.6 


1.93 


[Arm] 7136+.. 


23.5 


23.5 


1.00 


23.5 


1.00 


[Arm] 8.99 /im 




8.6 




8.1 


- 


[Ariv] 4711 1 


8.6 


1.5 


0.76 


8.2 


1.53 


Hei4713 1 




5.0 




5.0 




[Ariv] 4740 


4.5 


1.2 


0.26 


6.2 


1.39 


[Fell] 5.34 yum 




0.1 




10.8 




[Fell] 26.0 yum 


34. 


0.0 


0.00 


3.4 


0.10 


[Fein] 4658 


4.5 


4.5 


1.00 


4.5 


1.00 


[Fein] 4986 


7.4 


5.8 


0.78 


7.0 


0.94 


[Fein] 22.9 jim. 




2.2 




3.2 




[Feiv] 4906 


3.0: 


1.6 


0.54 


3.1 


1.02 


[Fev] 4227 


1.8 


1.5 


0.84 


5.5 


3.10 



" In Col. 1, blends are indicated by braces. The observed intensity of a 
blend is attributed to the first line and an asterisk to the second line. 
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Table 4. Model 'predictions' for IZw 18NW: model intensities 
divided by observed intensities 



Line ident. 


Obs. 


NO- 


-Nl 


Ml 


M2 


M3 


M4 




2 


3 - 


-4 


5 


6 


7 


8 


eabs/e" 


0.37 


.20- 


-.43 


0.20 


0.34 


0.38 


0.43 


Civ 1549 


500: 


.14- 


-.41 


1.12 


0.65 


0.64 


1.08 


[Oi] 6300+ 


8.5 


.01- 


-.59 


1.19 


1.01 


0.95 


1.38 


[Oii] 7320+ 


6.3: 


1.5- 


1.24 


1.19 


1.18 


1.19 


1.15 


[Oiii] 4363 


65.9 


.73- 


-.82 


0.96 


0.96 


1.00 


0.94 


[Oiv] 25.9yu 


49.1 


.37- 


-.84 


1.08 


0.97 


0.94 


0.98 


[Neil] 12.8yu 


9.0: 


.14- 


-.10 


0.14 


0.21 


0.22 


0.23 


[Neiii]15.5/i 


45.7 


1.31- 


-1.22 


1.09 


1.07 


1.02 


1.11 


[Sin] 34.8;u 


157. 


.03- 


-.10 


0.17 


0.14 


0.14 


0.18 


Sim] 1882+ 


270: 


.85- 


-.99 


1.12 


1.26 


1.34 


1.16 


[Sii] 6716 


22.5 


.30- 


-.54 


0.96 


0.78 


0.77 


1.10 


[Sii] 6731 


16.9 


.30- 


-.51 


0.96 


0.78 


0.77 


1.10 


[S II] 4068* 


3.7* 


.30- 


-.43 


0.73 


0.58 


0.58 


0.80 


[Siii] 9531+ 


114. 


1.15- 


-1.05 


0.95 


0.99 


1.00 


0.92 


[a III] 6312 


6.7 


.89- 


-.84 


0.77 


0.85 


0.87 


0.75 


[Siii] 18.7;/ 


28. 


1.15- 


-1.01 


0.94 


0.95 


0.94 


0.91 


[S III] 33.5yu 


120. 


.45- 


-.44 


0.39 


0.40 


0.40 


0.39 


[Siv] 10.5yU 


48. 


.87- 


1.50 


2.15 


1.93 


1.97 


2.15 


[Ariv] 4740 


4.5 


.26- 


-.73 


1.87 


1.39 


1.47 


1.76 


[Fell] 26.0yu 


34. 


.00- 


-.05 


0.10 


0.10 


0.10 


0.14 


[Feiii] 4986 


7.4 


.78- 


1.15 


0.90 


0.94 


0.89 


0.89 


[Feiv] 4906 


3.0: 


.54- 


-.70 


0.96 


1.02 


1.09 


1.02 


[Fev] 4227 


1.8: 


.84- 


-2.3 


2.1 


3.1 


3.2 


3.1 



"2abs/2- absolute values. 
''/14068 intensity not corrected for [Fe v]^4071 (Tabled. 



Table 5. Influence of shell parameters on selected line intensity 
predictions for models M2 and M4 



Param. / line 


M2 


M2b 




M2c 


M4b 


M4 


M4c 




2 


3 




4 


5 


6 


7 


fCOV 

J\ 


0.26 


0.26 




0.22 


0.29 


0.29 


0.22 


fcov 

Jl 


0.30 


0.60 




0.60 


0.30 


0.60 


0.60 


T„,(2) 


1.21 


0.62 




0.92 


2.17 


1.01 


1.57 




4.9 


4.7 




3.8 


4.6 


4.3 


3.0 


PJk/lO' 


3.4 


3.4 




3.1 


2.7 


2.7 


2.4 




21.7 


22.5 




24.6 


26.8 


26.8 


26.8 


O/HxlO^ 


168. 


165. 




162. 


178. 


173. 


171. 


[Oi] 6300+ 


1.01 


1.00 




0.92 


1.42 


1.37 


1.06 


[Oiii] 4363 


0.96 


0.97 




1.00 


0.92 


0.94 


0.95 


[Oiv] 25.9^ 


0.97 


0.91 




0.92 


1.05 


0.98 


0.99 


[Sii] 6716 


0.78 


0.81 




0.79 


1.05 


1.10 


0.84 


[S IV] I0.5fi 


1.93 


1.97 




2.05 


2.07 


2.15 


1.92 


[Ariv] 4740 


1.39 


1.48 




1.62 


1.55 


1.77 


1.81 


Table 6. Fine-structure line intensities" and X-rays 


Line* 


Obs. 


M2 


M2x 


M2x2 


H I colls.' 




no 


yes 


no 


yes 


no 


yes 


[Cii] 158. 




2.1 


2.1 


124 


119 


224 


204 


[Oi] 63.2 




4.7 


4.6 


274 


210 


624 


450 


[Neil] 12.8 


9 


1.9 


3.6 


1.9 


5.5 


2.0 


16 


[Sin] 34.8 


157'' 


21 


22 


37 


80 


55 


190 


[Am] 7.0 




1.2 


1.3 


1.4 


2.2 


2.2 


7.9 


[Fell] 26.0 


34 


3.4 


4.8 


5.8 


38 


9.3 


108 


[Fell] 35.3 




0.8 


1.4 


1.1 


9.5 


1.7 


33 



" In units HyS = 1000. 
* Wavelengths in yum. 

Inelastic collisions of [Sin], [Feii], [Neii] and [Am] with Hi omitted 
('no') or included ('yes'). 

'' The re-calibrated [Sin] flux is 68 after Sect. [6A2l 



gas distribution and photoionization by stars in IZw 18 

Table 7. Abundances in IZwlSNW compared to solar and 
Galactic halo stars 







X/O 


X/O 


[X/H] 


[X/Fe] 


[X/Fe] 


El. 


M" 


M 


IT* 


M 


M 


H' 


C 


6.55 


-0.66 


-0.77 


-1.84 


-0.02 


-0.05+0.15 


N 


5.60 


-1.61 


-1.56 


-2.18 


-0.36 


-0.40+0.30 





7.21 






-1.45 


0.37 


0.52+0.15 


Ne 


6.39 


-0.82 


-0.80 


-1.45 


0.37 




Si 


5.89 


-1.32 


-1.46 


-1.62 


0.20 


0.21+0.08 


S 


5.57 


-1.64 


-1.55 


-1.57 


0.25 


0.21+0.07 


Ar 


4.97 


-2.24 


-2.16 


-1.55 


0.27 




Fe 


5.63 


-1.58 


-1.45 


-1.82 







" 12 + log(X/H) by number in present model M. 

* Empirical abundance: Izotov & Thuan (1998, 1999) with 12 + 

log(0/H) = 7.16; Fe: Thuan & Izotov (2005). 

" Halo stars [Fe/H] ~ -1.8: Fabbian et al. (2006), Garcia Perez et al. 
(2006), Nissen et al. (2007), Carbon et al. (1987). 
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Table D.l. Effective collision sttengths for [O m] 



T* /I n4 IT. 


0.5 


1.0 


2.U 


3.0 


Reference:" 










1 T-\\ 

- D) 






Sea58 


- 




1 


59 


- 


- 




SSS69 


- 




2 


39 


- 


- 




ENS69 


1 


85 


2 


50 


2.91 


2 


96 




2 


17 


2 










Men83 


2 


02 


2 


17 


2.39 






Ag83 


2 


035 


2 


184 


lAQA 


2 


511 


BLS89 


2 


10 


2 


29 


2.51 


2 


60 


Ag93 


2 


039 


2 


191 


2.414 


2 


519 


LB 94 


2 


1268 


2 


2892 


2.5174 


2 


6190 


AgK99* 


2 


0385 


2 


1906 


2.4147 


2 


5191 


LB94/AgK99^ 


1 


0435 


1 


0450 


1.0425 


1 


0397 










n(3p 


- S) 






Sea58 


- 







220 


- 


- 




SSS69 


- 







335 


- 


- 




ENS69 





255 





298 


0.331 





339 


ES74 





276 





325 








Men83 





248 





276 


0.314 






Ag83 





2521 





2793 


0.3162 





3315 


BLS89 





260 





287 


0.318 





331 


Ag93 





2732 





2885 


0.3221 





3404 


LB94 





2720 





2925 


0.3290 





3466 


AgK99'' 





2732 





2885 


0.3221 





3404 


LB94/A2K99^ 





9956 


1 


0139 


1.0214 


1 


0182 












- S) 






Sea58 


- 







640 


- 


- 




SSS69 









310 


- 






ENS69 


n 
u 




u 


J / o 


0.555 






ES74 





807 





856 


0.752 






Men83 





516 





617 


0.634 






Ag83 





5463 





6468 


0.6670 





6524 


BLS89 





59 





677 


0.664 





634 


Ag93 





4312 





5227 


0.5769 





5812 


LB94 





4942 





5815 


0.6105 





6044 


AgK99* 





4312 





5227 


0.5769 





5812 


LB94/AgK99<^ 


1 


1461 


1 


1125 


1.0582 


1 


0399 



" Refs: Sea58: Seaton(1958); SSS69: Saraph et al. (1969); ENS69: 
Eissner et al. (1969); ES74: Eissner & Seaton (1974); Men83: Mendoza 
(1983); Ag83: Aggarwal (1983); BLS89: Burke et al. (1989); Ag93: 
Aggarwal (1993); LB94: Lennon & Burke (1994); AgK99: Aggarwal 
& Keenan (1999). 
* Results from Aggarwal (1993). 
" Collision sttength ratio. 



